Cascades and Collapses, Great Walls and Forbidden Cities: 
Infinite Towers of Metastable Vacua in Supersymmetric Field Theories 
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In this paper, we present a series of supersymmetric models exhibiting an entirely new vacuum 
structure: towers of metastable vacua with higher and higher energies. As the number of vacua 
grows towards infinity, the energy of the highest vacuum remains fixed while the energy of the true 
ground state tends towards zero. We study the instanton-induced tunneling dynamics associated 
with such vacuum towers, and find that many distinct decay patterns along the tower are possible: 
these include not only regions of vacua experiencing direct collapses and/or tumbling cascades, 
but also other regions of vacua whose stability is protected by "great walls" as well as regions of 
vacua populating "forbidden cities" into which tunnelling cannot occur. We also discuss possible 
applications of this setup for the cosmological-constant problem, for studies of the string landscape, 
for supersymmetry breaking, and for Z' phenomenology. Finally, we point out that a limiting case 
of our setup yields theories with yet another new vacuum structure: infinite numbers of degenerate 
vacua. As a result, the true ground states of such theories are Bloch waves, with energy eigenvalues 
approximating a continuum and giving rise to a vacuum "band" structure. 



I. INTRODUCTION 

The vacuum structure of any physical theory plays 
a significant and often crucial role in determining the 
physical properties of that theory. Indeed, critical issues 
such as the presence or absence of spontaneous symmetry 
breaking often depend entirely on the vacuum structure 
of the theory in question. 

Likewise, it may happen that a given model contains 
not only a true ground state, but also a metastable vac- 
uum state above it. Such models are also of considerable 
interest, for even when the true ground state preserves 
the apparent symmetries of a model, the physical proper- 
ties associated with the metastable vacua can often differ 
markedly from those of the ground state. In such situ- 
ations, the resulting phenomenology of the model might 
be determined by the properties of a metastable vacuum 
rather than by those of the true ground state. 

In recent years, models containing metastable vacua 
have captured considered attention. This is true for a va- 
riety of reasons. For example, metastable vacua can serve 
as a tool for breaking supersymmetry P, Q in the con- 
text of certain supersymmetric non-Abelian gauge theo- 
ries which are otherwise known to contain supersymmet- 
ric ground states. In addition, theories with large num- 
bers of vacua have also been exploited in various ways as 
a means of addressing the cosmological-constant prob- 
lem 0, 0) H @1 and obtaining de Sitter vacua in string 
compactifications (3]. Furthermore, the possibility of 
phase transitions in theories with multiple (meta)stable 
vacua leads to a number of implications for cosmology. 

Such ideas provide ample motivation to investigate 
whether there might exist relatively simple field theo- 
ries which give rise to additional, heretofore-unexplored 
vacuum structures. If so, such structures could poten- 
tially provide new ways of addressing a variety of un- 
solved questions about the universe we inhabit. 

In this and a subsequent companion paper Q, we will 



demonstrate that two new non-trivial vacuum structures 
are possible in relatively simple supersymmetric field the- 
ories. Moreover, as we shall see, the models which give 
rise to these non-trivial vacuum structures are not es- 
oteric; they are, in fact, simple generalizations of C/(l) 
quiver gauge theories. 

• First, we shall demonstrate through an explicit con- 
struction that certain supersymmetric field theo- 
ries can give rise to large (and even infinite) tow- 
ers of metastable vacua with higher and higher en- 
ergies. The emergence and analysis of this vac- 
uum structure will be the primary focus of the 
present work. As we shall see, as the number of 
vacua grows towards infinity in such models, the 
energy of the highest vacuum remains fixed while 
the energy of the true ground state tends towards 
zero. We shall study the instanton-induced tunnel- 
ing dynamics associated with such vacuum towers, 
and find that many distinct decay patterns along 
the tower are possible: these include not only re- 
gions of vacua experiencing direct collapses and/or 
tumbling cascades, but also other regions of vacua 
whose stability is protected by "great walls" as well 
as regions of vacua populating "forbidden cities" 
into which tunnelling cannot occur. Furthermore, 
as we shall see, these vacua are phenomenologically 
distinct from one another in terms of their mass 
spectra and effective interactions. 

• Second, we shall also show that there exists a lim- 
iting case of the above construction in which all of 
these infinite metastable vacua become degenerate, 
and in which there emerges a shift symmetry re- 
lating one vacuum to the next. As a result, the 
true ground states of such theories are nothing but 
Bloch waves across these degenerate ground states, 
with energy eigenvalues approximating a contin- 
uum and giving rise to a vacuum "band" structure. 
In this paper, we will merely sketch how such a 



2 



vacuum structure emerges; the complete analysis 
of such a structure will be the subject of a subse- 
quent companion paper [1]. 

This paper is organized as follows. In Sect. II, we 
present the framework on which our model is based. As 
we shall see, our model is essentially nothing more than 
an Abelian quiver gauge theory, expanded to allow ki- 
netic mixing between the various U{1) factors. We shall 
then proceed to discuss the corresponding vacuum struc- 
ture which emerges from this framework, including all 
stable vacua and all saddle-point barriers between them. 
We shall also discuss radiative corrections to this vacuum 
structure, and demonstrate that these corrections can be 
kept under control. In Sect. Ill, we shall then discuss the 
decay dynamics along these metastable vacuum towers, 
and examine the different sorts of instanton-induced tun- 
neling decay patterns which are possible. In Sect. IV, wc 
then analyze the particle spectra in each vacuum of the 
tower, and demonstrate how these spectra evolve as our 
system tumbles down the vacuum tower. In Sect. V we 
shift gears briefly, and consider the limiting case of our 
scenario in which our infinite towers of metastable vacua 
become an infinite series of degenerate ground states. 
Thus, in this limit, the true ground states of such theories 
are Blocli waves. Finally, in Sect. VI, we enumerate the 
potential physical applications of our vacuum towers, in- 
cluding possible new ideas for the cosmological-constant 
problem, for studies of the string landscape, and for Z' 
phenomenology. 

We emphasize that our primary goal in both papers is 
the demonstration that such non-trivial vacuum struc- 
tures can emerge in relatively simple supersymmetric 
field theories. Although there exist numerous implica- 
tions and applications of these ideas (some of which will 
be discussed in Sect. VI), our primary goal in these pa- 
pers will be the study of the emergence and properties of 
these vacuum structures themselves. 



II. THE GENERAL FRAMEWORK 

We begin by presenting our series of supersymmetric 
models which give rise to infinite towers of metastable 
vacua. Specifically, for each N > 1, we shall present a 
model which contains not only a stable vacuum ground 
state but also a tower of iV — 2 metastable vacua above it. 
Our model consists of N different t/(l) gauge group fac- 
tors, denoted U{l)a (a = 1, ■■■,N), as well as A^-|-l differ- 
ent chiral superfields, denoted $i {i = 1, A^-l-l). These 
superfields carry the U{1) charge assignments shown in 
Table [H and follow the well-known orbifolded "moose" 
pattern wherein each field with 2 < i < N simultane- 
ously carries both a positive unit charge under U{l)i-i 
and a negative unit charge under U{l)i. By contrast, 
the fields <I>i and ^n+i sit at the orbifold endpoints 
of the moose, and are charged only under the corre- 
sponding endpoint gauge groups respectively. We shall 
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TABLE I: The field content and charge assignments for the 
chiral superfields in the models under consideration. 

assume that each U{l)a gauge field has a correspond- 
ing gauge coupling g^, and for simplicity we shall fur- 
ther assume that ga = g for all a. Note that the only 
non-vanishing gauge anomalies inherent in this charge 
configuration are mixed C/(l) anomalies proportional to 
J^i'^^ QaiQbij which can be canceled by the variant of 
the Green-Schwarz mechanism Q discussed in Ref. [loj . 

To this core model we then add three critical ingre- 
dients, each of which is vital for the emergence of our 
metastable vacuum towers. First, given the field content 
of each model, we see that the most general superpoten- 
tial that can be formed in each case consists of a single 
Wilson-line operator 

iV+l 

= A [| . (1) 

We shall therefore assume that this operator is turned on 
for each value of TV. Note that the coupling A has mass 
dimension 2 — N and is therefore non-renormalizable for 
all TV > 2. 

Our second and third ingredients both exploit the 
Abelian nature of our gauge groups. The second ingre- 
dient is to introduce non-zero Fayet-Iliopoulos terms i^i 
and for the "endpoint" gauge groups C/(l)i and U{1)n 
respectively. While all of our C/(l)a gauge groups could 
in principle have corresponding non-zero Fayet-Iliopoulos 
terms , we shall see that turning on only and will 
be sufficient for our purposes. Indeed, we shall prune our 
model further by taking = = ^. 

Finally, our third ingredient is a simple one: kinetic 
mixing 11]. It is well-known that the field strength ten- 
sor F'^'^ for an Abelian gauge group is gauge invariant by 
itself. Thus, in a theory involving multiple U{1) gauge 
groups, nothing forbids terms proportional to F^^Fi,,^i, 
from appearing as kinetic terms in the Lagrangian, where 
a ^ h. Similarly, in a supersymmetric model, mixing 
between the field-strength superfields is permitted, 
whereupon the gauge-kinetic part of the Lagrangian may 
take the generic form [l2| 

C 3 ^ [d^e Wac. Xab (2) 
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with a general (symmetric) kinetic-mixing matrix Xat- 
I 



^ah — 



1 

-Xl3 



-Xl2 
1 

-X23 



-Xl3 
-X23 
1 



-X2Af -X3Af 



-XiJv\ 

-X2Ar 
-X3Ar 

1 / 



(3) 



As long as Xab is non-singular, with positive real eigen- 
values, there exists a matrix Mab which transforms the 
U(l) gauge groups into a basis in which their gauge- 
kinetic terms are diagonal and canonically normalized. 
Specifically, one can write WjXabWb = (Wa)^(W^o) 
where Wa = MabWb- In general, such a matrix M takes 
the form M = SO where 5 is a diagonal rescaling matrix 
whose entries are the square roots of the eigenvalues of 
X and where O is an orthogonal rotation matrix diago- 
nalizing X. After this diagonalization process, the new 
U{l)a charge assignments Qai for our fields $i and the 
new Fayet-Iliopoulos parameters for our gauge groups 
U{l)a are given in terms of the quantities Qai and h^ 
the original basis: 



Qa^ = [{M-yUQb 



(4) 



In this vein, it is important to note that the matrix M 
corresponding to each kinetic-mixing matrix X is not 
unique. Any matrix of the form M' ~ VM, where V 
is an orthogonal matrix, also yields the correct normal- 
ization for the gauge-kinetic terms. Different choices for 
V correspond to different orthogonal choices for the fi- 
nal basis of C/(l)'s. Ultimately, the physics is insensitive 
to which basis is chosen. By contrast, the rescaling ma- 
trix S is unique, and it is this matrix which carries the 
physical effects of kinetic mixing. 

In general, any of the Xab parameters in Eq. ([3]) may 
be non-zero. However, it will be sufficient for our pur- 
poses to restrict our attention to the case in which only 
"nearest- neighbor" [/(l)'s experience mixing. Thus we 
shall assume that Xab ^ if and only if 6 = a -I- 1. For 
simplicity, we shall further assume that all non-zero Xab 
are equal, so that Xab — X^a+i,b- While more general 
kinetic-mixing parameters may be chosen, we shall see 
that these simplifications enable us to expose the exis- 
tence of our metastable vacuum towers most directly. 

Needless to say, it would have been possible to con- 
struct our models entirely without kinetic mixing by pos- 
tulating highly non-trivial choices for Qai and right 
from the beginning. However, we have found that it is 
easier to begin with the simpler assignments Qai shown 
in Table U and the Fayet-Iliopoulos terms = ^jv = ^ 
described above, and to bundle all of the remaining com- 
plexities in terms of a single kinetic- mixing parameter 
X- 

Not all values of the parameter x lead to self-consistent 
theories, however; we must also ensure that the kinetic- 
mixing matrix Xab in Eq. ([3]) is invertible with positive 




FIG. 1: The maximum allowed value of |x| 
tion of A*'. As TV — ^ oo, we see that |x|n 
approaches 1/2 from above. 



plotted as a func- 
ax asymptotically 



(real) eigenvalues. For N — 2, we find that this requires 
Ixl < 1, while for = 3, 4, and 5 this requires |x| < 
1/V2, Ixl < 2/(1 + V5), and |x| < 1/^3 respectively. 
The behavior of the maximum allowed value of |x| as a 
function of N is shown in Fig. [TJ 

For arbitrary N, we see from Fig.[T]that the maximum 
allowed value of |x| always exceeds 1/2, and asymptoti- 
cally approaches 1/2 as A^ oo. Moreover, we find that 
negative values of x do not lead to the metastable vac- 
uum towers which are our main interest in this paper. 
As a result, we shall simplify matters by restricting our 
attention to the range 



< X < 1/2 



(5) 



(Indeed, only in Sect. V shall we consider the x = 1/2 
limit.) Likewise, our orbifold moose structure for any 
length A^ possesses a reflection symmetry under which 
the combined transformations S^a —^N+i-a, 9a ~* 



9N+1- 



6 



-^2, and <bj <& 



N-j+2 



leaves the physics 



invariant. This means that the scalar potential in a the- 
ory of given N with ^ < will be identical to that with 
^ > 0, save that the role played by $i in the former is 
played by ^n+i in the latter, and so forth. As a result, 
we will restrict our attention to situations with 



C > 0. 



(6) 



Finally, our model also has a reflection symmetry un- 
der A ^ —A, as a result of which we can further re- 
strict to A > 0. However, for each N , we shall flnd that 
there is actually a minimum positive value A^ which is 
needed in order for our entire tower of A^ — 1 vacua to be 
(meta)stable. The derivation and interpretation of this 
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critical value will be discussed further below. We 
shall therefore actually restrict to the range 

A > (7) 

in much of what follows. 

Thus, to summarize, our models are defined in terms 
of N different U{l)a gauge groups and iV + 1 different 
chiral superfields $i , with charges indicated in Table [J 
This structure can be indicated pictorially through the 
moose diagram in Fig. [21 which shows not only the U{l)a 
gauge groups but also the <i>i fields which provide nearest- 
neighbor "links" between them as well as the x parame- 
ter which governs their universal nearest-neighbor kinetic 
mixing. For each value of TV > 2, our models are there- 
fore governed by four continuous parameters, namely g, 
X, ^, and A, subject to the bounds in Eqs. (O, ® and 
([7]). Note that a similar model, but with x = 0, was 
considered in Ref. ,13]. 

FIG. 2: The "moose" diagram for our series of models. The 
A'^ sites [each representing a (7(1) gauge group] are connected 
by links corresponding to the chiral superfields <l>i , and experi- 
ence nearest-neighbor kinetic mixing governed by a universal 
parameter x- 

Our main interest in this paper is in the vacuum struc- 
ture of these models. This in turn is governed by their 
corresponding scalar potentials. In general, the scalar po- 
tential for such a supersymmetric gauge theory coupled 
to matter includes both I?-tcrm and i^-term contribu- 
tions and can be written in the form 

N N+1 

V = lY.g^Dl+J2\F,\\ (8) 

a—l i—1 

where each gauge-group factor has a common coupling g 
and where 

For any choice of parameters {5, x, ^, A}, the extrema 
of the scalar potential can then be obtained by solving 
the TV -|- 1 coupled simultaneous equations 

dV 

— = (* = l,...,iV+l) . (10) 

d(t>i 

However, a solution is a local minimum only if the eigen- 
values of the 2(iV+ 1) X 2(iV+ 1) mass matrix 



(11) 



are all non-negative and the number of zero eigenvalues is 
precisely equal to the number of Goldstone bosons eaten 
by the massive gauge fields. (Indeed, additional zeroes 
would indicate the presence of classical fiat directions.) 
In what follows, however, we will use the term "vacuum" 
loosely to refer to any extremum of the potential and em- 
ploy adjectives such as "stable" and "unstable" to distin- 
guish the eigenvalues of the mass matrix. Of course, a 
"metastable" vacuum exists only when two or more vacua 
exist and are stable according to the above definitions; 
all but the vacuum with lowest energy are considered 
metastable. 

Note that the scalar potential V will in general be a 
function of only the absolute squares of fields. As a result, 
we can take all non-zero vacuum expectation values Vi = 
{(pi) to be real and positive without loss of generality. 

Also note that while the F-terms are insensitive to the 
kinetic mixing, the Z?-terms in Eq. ([9]) are calculated in 
terms of the charges Qai and Fayet-Iliopoulos coefficients 

in the orthonormalized basis given in Eq. ([3]) . It is ap- 
parent from the form of Eq. ^ that while Da depends on 
the choice of O within the matrix M, the scalar poten- 
tial as a whole is insensitive to this choice. By contrast, 
the rescaling matrix S within M is physical, modifying 
the Z?-terms in a non-trivial way. It is in this manner 
that the effects of kinetic mixing are felt in the vacuum 
structure of the theory. 

Finally, we emphasize that we shall deem an extremum 
of the scalar potential V to be (meta)stable only if there 
are neither fiat directions nor negative eigenvalues in the 
mass matrix. These two restrictions are rather severe, 
since most models tend to give extrema which have either 
tachyonic modes or fiat directions. 

In the rest of this paper, we shall simplify our anal- 
ysis by scaling out the gauge coupling g in the manner 
discussed in Ref. [13] ■ Specifically, we shall define the 
rescaled quantities ^' = 5^, A' = A/.g^/^, and $■ = y^$i, 
holding all other quantities fixed. We shall then eliminate 
explicit dependence on ^' by further rescaling all dimen- 
sionful quantities by appropriate powers of ^' in order to 
render them dimensionless. Specifically, we shall define 
A" = A7(e')'"^/', = $,/VF, and V" = V/{^')^. 
In practical terms, the net effect of these two rescal- 
ings is that we simply rewrite all of our original expres- 
sions in terms of the new variables A" = £,^^'^~^\/ g, 
V" = V/{g£,f, and = and then drop the 

double primes. Thus, for each N >2, our models can be 
analyzed purely in terms of a single kinetic- mixing pa- 
rameter X and the rescaled (dimensionless) Wilson-line 
coefficient A defined above; the resulting vacuum energies 
V and field VEV's Vi = will then be dimensionless 
as well. Finally, we shall adopt a notation (first intro- 
duced in Ref. [13| ) wherein we describe a particular field 
configuration of VEV's Vi as belonging to a class denoted 
{p, q, ...} if the only non-zero VEV's for the vacuum so- 
lutions in this class are Vp, Vq, and so forth. For example, 
{1,2,4} will refer to a vacuum configuration in which wi, 
V2, and U4 are non-zero, with all other VEV's vanishing. 
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Our claim, then, is that for each N, our model gives 
rise to a tower consisting of A'^ — 1 vacuum solutions. 
Specifically, for each N, we claim that the corresponding 
model will give rise to a true stable ground state along 
with a series of — 2 metastable ground states with 
higher and higher vacuum energies. 



{123} 

as X- 



{1} 

' 00 



N=3 



n=1 
{12} 



A. Example: iV = 3 

We shall begin the analysis of our models by focus- 
ing on a simple example: the = 3 special case, which 
consists of three U{1) gauge groups and four chiral su- 
perfields. The scalar potential in this case is given by 
V = Vd + Vf, where 

Vd = l{\<Pl\'-\^2f-\h\' + \M')' 

.|0,|2„(V2-1)|02P 

(V2-l)|(/)2|2 + |^4|^ 

2~ |</)ip + (\/2+l)|02|^ 



8(1 - V2x) 



8(1 + V2x) 



and 



Vf = |A| 



4 

1=1 



(V2+l)|(/)2p + |</'4 



, (12) 



(13) 



In these expressions, of course, 0i denotes the (complex) 
scalar component of the chiral supermultiplct $i. 

Given the scalar potential, it is then a straightforward 
matter to calculate the vacuum structure of this poten- 
tial. Our results are as follows. Defining 



1 



(14) 



we find that there are two (meta)stable vacua in this 
model for all A > A3. The first vacuum state (which we 
shall call the n — 1 vacuum) has energy Vi = 1/2 and 
corresponds to the solution with 



n = 1 



vl = 1 + X 



V2=X 



vl^vl^O . (15) 



By contrast, the second vacuum state (which we shall call 
the n = 2 vacuum) has vacuum energy V2 = — x)^^ 
and corresponds to the solution with 



n = 2 : 



X 



2x 



vl = 



(16) 

We thus see that the n ~ 1 vacuum is of {12}-type, 
while the n = 2 vacuum is of {13}-type. As emphasized 



{13} 



FIG. 3: A sketch of the vacuum structure of the A'^ = 3 
model. For A > A3, the corresponding scalar potential gives 
rise to two distinct minima: a {13} vacuum which serves as 
the true ground state of the theory, and a {12} vacuum which 
serves as an additional, metastable vacuum. These two min- 
ima are separated by a saddle-point {123} extremum which 
reduces to the {1} extremum in the formal A — » 00 limit. Note 
that this sketch is actually a two-dimensional representation 
of potential energy contours in a three-dimensional field space 
parametrized by {vl,V2,v^}. 



above, both of these solutions are stable (without any flat 
or tachyonic directions) for all A > A3; however, since 
Vi > V2 for all X < 1/2, we see that the n — 2 vacuum 
is the true ground state in this theory, while the n ^ 1 
vacuum is only metastable. This vacuum configuration 
is sketched in Fig. [31 

Note that for A > A3, the n — 1 and n = 2 vacua 
are separated by a potential barrier whose lowest point 
is a {123} saddle-point extremum of the scalar poten- 
tial. Unlike the field-space solutions for the n = 1 and 
n ~ 2 vacua, which are A-independent, the field-space 
solution for this saddle point depends quite strongly on 
A. This is shown in Fig. 31 where the explicit solutions 
for {vi,V2,vl} are plotted for A > A3. 

It will be useful to understand how this vacuum struc- 
ture deforms as a function of A. Naively, one might sus- 
pect that taking X —^ 00 would cause the height of this 
saddle-point barrier to diverge. However, we see from 
Fig. m that this is not the case: the scalar potential V 
and all of the field VEV's vf quickly reach finite asymp- 
totes. Indeed, in the formal A ^ 00 limit, we see that 



"3 



0, whereupon our {123} saddle-point solution 



reduces to the {1} solution given by 



1 



2 2 2 n 



(17) 



with 1/12 = i(l-x')-^ 

The above results are valid for all A > A3. However, 
it will also be important for us to understand what hap- 
pens as we reduce the value of A below A3. Since A sets 
the scale for the barrier height between the two vacua in 
Fig. 121 reducing A has the effect of reducing the bar- 
rier height between the two vacua. This in turn will 
destabilize our metastable vacuum. Specifically, we see 
from Fig. [H that A3 = l/-\/x(l + x) is nothing but the 
critical value of A at which the barrier height becomes 



X 



X 



0,525 I 

N = 3 




FIG. 4: The {123} saddle-point solution for A'^ = 3 and x ~ 1/5, plotted as a function of A > A3 = 5/VG ~ 2.04. Here we show 
the non-zero field VEV's (upper left plot), \v2\'^ (upper right plot), and \v3\^ (lower left plot), as well as the corresponding 
scalar potential V (lower right plot). While each of these quantities varies with A, they quickly reach formal asymptotes as 
A 00. 



equal to Vi. At this critical value, the {123} saddle- 
point solution shown in Fig. |3] actually merges with the 
metastable n = 1 vacuum solution (which is of {12}- 
type) and thereby destabilizes it. Thus, as we take A 
below A3, we lose the n = 1 vacuum. In order to empha- 
size that this is the critical A-value at which the n = 1 
vacuum is destabilized, we shall also refer to A3 as A3 . 

Taking A still lower, we ultimately reach a second crit- 
ical value A3 2 = ^^2/(2 — x) at which even the n — 2 
vacuum becomes destabilized. In this case, a new {123} 
solution becomes stable and serves as the ground state 
of the theory for all A < A3 2- 

These two critical values of A for the = 3 case are 
plotted in Fig. [5] as a function of the kinetic-mixing pa- 
rameter X- Note, in particular, that A3 i diverges as 
X ^ 0. This shows that the stability of the n — 1 



metastable vacuum solution relies not only on having a 
sufficiently large value of A, but also on the existence of 
non-zero kinetic mixing. Note also that 

< A3 2 < A* i for all < x < 1/2 . (18) 

This indicates that as reduce A below A3 = -^31, our 
= 3 metastable "tower" destabilizes from the top 
down, with the n = 1 metastable vacuum destabilizing 
before the n = 2 vacuum. It is for this reason that we 
can associate A3 (the critical A-value for stability for the 
entire tower) with A3 i (the critical A-value for stability 
of the highest vacuum) . 

We see, then, that our N = 3 model gives rise to 
two vacuum solutions (one stable ground state and one 
metastable state above it) for all A > AJ = A3 . The 
solutions for these vacua are A-independent, and depend 
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0.5 



FIG. 5: Critical A-values for the = 3 model, plotted as 
functions of X- Here A3 — A34 — l/y^xix + 1) a-nd A3^2 ~ 
1^/2/(2 — x)- For A > Aa i, both the n — 1 and n = 2 vacua 
are stable, while for AJj ^ < A < A3 1, the n = 1 vacuum 
is destabilized and only the n = 2 vacuum is stable. For 
A < A3_2, even the n — 2 vacuum is destabilized; in this range 
a new {123} solution becomes stable and serves as the ground 
state of the theory. 



only on x- However, the barrier height between these 
two vacua (and hence the stabihty of the metastable vac- 
uum) depends intimately on both x and A, and formally 
reaches an asymptote as A ^ cxd. 



B. Example: A*^ = 4 

Having discussed the vacuum structure of the N = 
3 model, we now turn to the N = 4 model. In this 
case, there are four C/(l) gauge groups and five chiral 
superfields. Defining 



a: 



(19) 



we find that the vacuum structure now consists of three 
vacuum solutions, each without tachyonic or flat direc- 
tions, for all A > A4. The n = I vacuum has energy 
Vi = 1/2, just as in the iV = 3 case, and corresponds to 
the solution with 




X 



(20) 



while the n ~ 2 vacuum has energy V = j{l — x)~^ i 
again just as in the = 3 case, and now corresponds to 
the solution with 

2-X 

' "2 — 2-2x ' 

(21) 



2-2x ' 



The solutions in Eqs. ([201) and ([21]) are clearly the iV = 4 
generalizations of the corresponding A'^ = 3 solutions in 
Eqs. (fT5| and (fT6|) respectively. However, the crucial new 
feature for the = 4 case is the appearance of an addi- 
tional vacuum state of even lower energy. This vacuum 
state, which we shall refer to as the n = 3 vacuum, has 
energy V3 = ^(3 — 4%)"^ and corresponds to the solution 
with 



,,2 _ 3(l-x) 



V. 



2 _ 2-3x 
3~4x 



,,2 - 1-X 
-'S - 3-4x ' 

' = «2 ^ . 



(22) 



This vacuum structure is sketched in Fig. [G] Note that 
just as in the = 3 case, these vacuum solutions are all 
A- independent. 




{123} . 



{1234} ^{13} 
as X-^ 00 



{134} 



FIG. 6: A sketch of the vacuum structure of the A'' = 4 
model. For A > A4, the corresponding scalar potential gives 
rise to three distinct minima: a {134} vacuum which serves 
as the true ground state of the theory, and two additional 
metastable vacua of types {123} and {124} above it. These 
three vacua are separated by three different saddle-point ex- 
trema which are each different solutions of {1234}-type; in 
the formal A 00 limit, these reduce to the different two- 
VEV solutions shown. Note that this sketch is actually a two- 
dimensional representation of potential energy contours in a 
four-dimensional field space parametrized by {iii , nf, , «!}. 

For A > A4, the 71 = 1, n = 2, and n — 'i vacua are sep- 
arated by potential barriers whose lowest points are all 
{1234} saddle-point extrema of the scalar potential. Un- 
like the field-space solutions for the (meta)stable vacua, 
the field-space solutions for these saddle points depend 
quite strongly on A. However, in the formal A — *■ 00 limit, 
these solutions all quickly reach finite asymptotes. For 
example, the asymptotic saddle-point solution between 
the n = 1 and n = 2 vacua is of {12}-type and is given 

by 



(1,2) 



i-x 



7)2 



3^ , 



1-X 



= = V? = 



(23) 



with energy V12 — ■^{1 — x'^) ^- This is the A'' = 4 ana- 
logue of the A^ = 3 asymptotic saddle-point solution in 
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Eq. p?)) . and continues to have the same energy Vu- 
However, in the = 4 case there are also additional 
saddle-point solutions which involve the new n — 3 vac- 
uum. Specifically, the asymptotic saddle-point solution 
which lies directly between the n = 1 and n = 3 vacua is 
of {13}-type and is given by 



9 _ 2(l-x) 



(1,3) 



x(i-x) 



(24) 



7)2 







with energy V13 = (1 — x)/(2 — 2^ — X^), while the asymp- 
totic saddle-point solution between the n = 2 and n = 3 
vacua is of {14}-type and is given by 



(2,3): 



2 _ 2(l-x) 
Jl - 2-2X-X' 



= vi = vi = 



„2 _ 1-X-X 



2-2X-X" 



(25) 



with energy V23 = i(2 - 2x - x')"'- 

This vacuum structure emerges for all A > A|. How- 
ever, just as in the = 3 case, we find that reducing A 
below A| tends to destabilize our vacuum tower. Specifi- 
cally, one finds that there are now three critical A-values, 



^^4,2, 



and A4 3 , at which the n = 1 , 



denoted A| 

and n = 3 vacua are respectively destabilized, 
three critical values are given by 



n = 2, 
These 



a: 



1 



4,1 



^4,2 = 



\* — 




3(2 -3x) 



(26) 



and all three are plotted in Fig. [7] as functions of x- We 
see from this figure that 



< At o < At 9 < At 



(27) 



for all < X < 1/2- This implies that just as in the 
iV = 3 case, our vacuum tower destabilizes from the top 
down as we reduce A below A| As a result of Eq. (P7)). 
we see that A| {i.e., the critical A-value for stabilizing 
the entire iV = 4 vacuum tower) is nothing but A4 ^ (the 
maximum of the individual values A| „ for stabilizing any 
of the individual vacua in the tower). It is this observa- 
tion which underlies the identification given in Eq. (jl9p . 
We also observe from Fig. [7] that non-zero kinetic mixing 
is also required in order for the stability of our metastable 
vacua. 

It is clear that the = 4 case is a direct generalization 
of the = 3 case. Once again, the vacuum solutions 
are A-independent, while the barrier heights (and hence 
the stability of these vacuum solutions) are A-dependent. 
Moreover, the energies associated with the n = 1 and 
n = 2 vacua, as well as the barrier between them, are 
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FIG. 7; Critical A-values for the A'^ = 4 model, plotted as 
functions of x- For A > A4 = all three of our vacua 

are stable, while for A4 2 < A < X41, the n = 1 vacuum is 
destabilized and for A2,3 < A < A2 2) both the n = 1 and 
n — 2 vacua are destabilized. Finally, for A < A4^3, even the 
n — 1 vacuum is destabilized; in this range a new {1234} 
solution becomes stable and serves as the ground state of the 
theory. 



unchanged in passing from the A^ = 3 case to the A^ = 4 
case. Indeed, the primary new feature in passing from 
the A' = 3 case to the A^ = 4 case is the emergence of a 
new vacuum solution, our so-called n = 3 vacuum, which 
"slides in" below the previous bottom of the tower and 
becomes the new ground state of the theory for all A > 
A| 3. As a result, our previous n — 2 ground state in the 
A^ = 3 theory now becomes the first-excited metastable 
state in the A^ = 4 theory, and our tower of metastable 
vacua has grown by one additional metastable vacuum. 
We stress again that all of these vacua are either strictly 
stable or strictly metastable. In particular, they contain 
neither tachyonic masses nor flat directions. 



C. Results for general A'^ 

We now turn to the case with general A^. As might be 
expected, the pattern we have seen in passing from the 
A^ = 3 case to the A^ = 4 case continues without major 
alteration. For general A^ and for A exceeding a critical 
value A^ , we find that our model has a vacuum structure 
consisting of a tower of A^ — 1 stable vacua: a single stable 
ground state, and N — 2 metastable vacua above it. As 
in the A^ = 3 and A^ = 4 cases discussed above, we shall 
number these vacua from the top down with an index n, 
so that the n = 1 vacuum sits at the top of the tower 
and the n = A^ — 1 vacuum (the true ground state of the 
theory) sits at the bottom. 
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We then find that the n-vacuum has energy 



where 



\xRr. 



l<n<N-l (28) 



(29) 



and corresponds to the solution with 



1 + l/R 





N+n- 



-1 - 1)1 Rn 



for i — 1 

for 2 < j < iV - n 
for j = TV - n + 1 
for TV - n + 2 < j < 
for j = N + 1 . 



(30) 



As evident from these VEV's, this is clearly a solution of 
{1, N — n, N — n + 2, N}-type. It is easy to verify 
that these results reduce to those already quoted for the 
N — i and iV = 4 special cases. Note that the vacuum 
energies Vn along the tower are independent of N , as 
expected; indeed, all that depends on N are the number 
of such vacua and their precise field VEV's. Also note 
that Vn-1 ^ as -/V ^ oo. 

As in the previous special cases with = 3 and = 4, 
any two vacua n and n' are separated by a saddle-point 
solution. In the formal A — *■ oo limit, we find that this 
barrier height asymptotes to the value 



1 



Rji' —n 



2x yRnRj^ 



- 1 



(31) 



where _R„ is defined in Eq. (j29p and where we have taken 
n' > n. Note that 14,™' is also independent of N. The cor- 
responding asymptotic saddle-point solutions are given 

by 



R^R„ 


R 



■j-N + n'-l' 
Rr^R„'-„- 







Rr,R„ 







for i — 1 

for 2 < j < - 

for j = A^ - n' + 1 

for A^ - n' + 2 < j < A^ - n 

for j = A^ — n -I- 1 

for A^ - n + 2 < j < Af 

for j = A^ + 1 . 

(32) 



These results are plotted in Fig. [5] for the A^ = 20 
model. Clearly, the vacuum structure of this model con- 
sists of a ground-state vacuum along with a tower of 18 
metastable vacua above it. In Fig. [SI we have shown the 
vacuum energies of these vacua, along with the asymp- 
totic energies of the saddle-point barriers which connect 



0.5 




5 10 15 

distance m field space 



FIG. 8: The vacuum structure of the Ai' = 20 model, plotted 
for X ~ 1/5. Each of the local minimum points corresponds to 
a (meta)stable vacuum state, while each local maximum point 
corresponds to the saddle-point configuration which directly 
connects the nearest-neighbor vacuum configurations on ei- 
ther side. For the purposes of this plot, we have then simply 
connected these points sequentially with straight lines. The 
vertical axis indicates the energies of these vacua or saddle 
points (the latter in their asymptotic limits), while the hor- 
izontal axis indicates the cumulative distances in field space 
along a trajectory which begins at the n — 1 vacuum and 
then proceeds along straight- line path segments to the (1,2) 
saddle point, then to the n = 2 vacuum, then to the (2,3) 
saddle point, and so forth. 



"nearest-neighbor" vacua, as functions of a cumulative 
distance in field space along a path that winds through 
each vacuum configuration and over each saddle point as 
it comes down the tower, vacuum by vacuum. In essence, 
then, this figure forms a linear "picture" of the tower. 

Despite the relative simplicity of Fig. [51 it is worth 
emphasizing that the geometry of the metastable vac- 
uum tower in the full (A'^ -I- l)-dimensional field space is 
rather non-trivial. For example, although the vacuum 
and saddle-point energies are plotted in Fig. [5| versus a 
cumulative, integrated distance as we wind our way down 
the vacuum tower, we could have just as easily defined 
the field-space distance associated with any vacuum or 
saddle point in terms of its straight-line distance directly 
back to a reference vacuum (such as the n = 1 vacuum 
at the top of the tower). The difference between these 
two different notions of distance is shown in Fig.[5|for the 
A^ = 20 case. Indeed, as evident from the actual solu- 
tions given Eq. (j30p , the vacua in our vacuum tower actu- 
ally lie along a "spiral" or "helix" in the A^-dimensional 
vn+1 ~ subspace of our full (A^ + l)-dimensional field 
space. 

It is also important to emphasize that Fig.[5|shows only 
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FIG. 9: The geometry of our metastable vacuum towers 
in field space, here illustrated for the A*' = 20 model with 
X ~ 1/5. For each n- vacuum (n = 1, 19) in the = 20 
vacuum tower, we have plotted two different notions of corre- 
sponding field-space distance. In curve (a), we have plotted 
the cumulative distance (as used in Fig. [S]) that accrues as 
we wind our way from the top of the tower down to the n- 
vacuum, passing through nearest-neighbor saddle points along 
the way. By contrast, in curve (b), we have plotted the direct 
straight-line distance between each n-vacuum and the refer- 
ence n = 1 vacuum at the top of the tower. The dramatic 
difference between these two curves is a reflection of the fact 
the vacua in our vacuum tower actually lie along a "spiral" 
or "helix" in the full 20-dimensional field space. 



those saddle-point barriers which exist between nearest- 
neighbor vacua along the vacuum tower. In actuality, 
however, there are saddle-point barriers which exist di- 
rectly between any two vacua (n,n'). As a result, it is 
possible to imagine descending through the vacuum tower 
taking "hops" with different values of An at each step. 
Fig. (TU] provides an illustration of the effects of different 
possiblities. 

Finally, we turn to the one remaining issue: the critical 
value „ at which the n}^ vacuum in the tower is desta- 
bilized. In general, for any N > 2 and l<n<A^— 1, it 



turns out that A 



\*2 



is given by 

m 



R 



N~2 



V 



r(nH-y) x{l + Rn) 



(33) 



where y = — 2%) = n/{Rn — 2) and where r(z) is 
the Euler F-function [for which F(z) ~ [z — 1)! when z G 
Z+J. It is straightforward to verify that these expressions 
reduce to the corresponding expressions for the A'^ = 3 
and N = A cases plotted in Figs. and [7] respectively. 
For example, using the result in Eq. (|33p . we find that 



A 



NA 



1 



X 



l-N/2 



(34) 



FIG. 10: The upper portions of the = 100 vacuum tower 
with X = 1/5, plotted for saddle-point "hops" of fixed magni- 
tudes (a) An = 1, (b) An = 3, and (c) An = 20 respectively. 
Note that the energy of a given vacuum is independent of how 
it is reached; for example, we see that the n = 4 vacuum has 
the same energy whether it is realized as the fourth minimum 
on the An = 1 curve or the second minimum on the An = 3 
curve. However, the corresponding field-space distances rela- 
tive to the top of the tower are generally smaller when larger 
hops are utilized when descending. 

In Fig. [nl we plot the values of A^ „ as functions of x 
for iV = 20 and 1 < n < 19. Unlike the simpler iV = 3 
and = 4 cases, however, we see that it is no longer 
true that A^ „ > AJ^r „/ for any n < n' . Instead, as N 
increases, we see that a complicated "crossing" pattern 
develops as a function of x- As a result of this crossing 
pattern, the value of n which results in the maximum 
value of A^ „ itself varies with x^ as shown in Fig. [T^ 
Nevertheless, for any value of x, we see that our entire 
vacuum tower will be stable if A > A^ , where 

(35) 



X* 



max 

l<n<Af- 



Thus, A^ corresponds to the upper "envelope" of the 
curves shown in Fig. 1111 

The crossing pattern shown in Fig. II II implies that our 
tower of metastable vacua will experience a non-trivial 
destabilization pattern as A is reduced from infinity. For 
any value of x in the range < x < 1/2, the first vacuum 
to be destabilized is indicated in Fig. [121 This serves as 
the initial destabilization point on the tower. Reducing 
A still further then induces a destabilization of the vacua 
immediately above and below this point, and further re- 
ductions in A result in a destabilization "wave" which si- 
multaneously runs both up and down the vacuum tower 
from this initial point until ultimately all vacua are desta- 
bilized. 

Finally, the dependence of A^ on N is shown in Fig.[T51 
We see that A^ generally grows rather quickly with N . 
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FIG. 11: Values of AJV,n/-^jvii plotted as functions of x for 
A'^ = 20 and 1 < n < 19. The curves corresponding to 
n = 1, 2, 4, 8, 19 are highlighted and labeled; the rest are not 
highlighted but proceed in sequential order relative to those 
which are. For x ~ 0.05, we see that the n = 2 vacuum has 
the highest value of „, while it is the n = 4 vacuum which 
has this property for x ~ 0.45. The upper "envelope" of all of 
these curves corresponds to the critical stability value for 
the vacuum tower as a whole . Note that all of these curves 
converge to A^^ „ = ^2^-^ j?, at x = 1/2. 
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FIG. 13: The critical values A^, plotted as functions of x for 
(a) N = 10, (b) iV = 30, (c) TV = 50, and (d) N = 100. We 
observe that A^ generally grows with TV, diverging as x ^ 
and asymptoting to ^J2'^~^ /Ji as x — > 1/2. We stress that the 
A-values plotted here are the rescaled, dimensionless versions 
of our original Wilson- line coefficients. As such, there is no 
conflict with the perturbativity of our model. 




FIG. 12: The n-index of the vacuum with the largest corre- 
sponding critical value AJf „, plotted as a function of x- This 
index specifies which vacuum in the tower is the first to desta- 
bilize as A is lowered from infinity. In general, for large TV, this 
index shifts from n = 1 to n = [TV/2] (where [x] is the greatest 
integer < a;) as x increases within the range < x < 1/2- 



Moreover, as already anticipated from the TV = 3 and 
TV = 4 special cases, A^ always diverges as x ~^ 
and asymptotes to A^^ ^2^^V3 as x ^ 1/2- The 
fact that A^ diverges as x ~* for every TV indicates 
that kinetic mixing plays a critical role in keeping our 
metastable vacuum tower stable. 

We also stress that the A-values plotted in Fig. [T3] 
are the rescaled, dimensionless versions of our original 
Wilson-line coefficients. Indeed, they are only rescaled 
variables, not to be confused with our primordial La- 
grangian couplings; indeed, as discussed earlier, this 
rescaling absorbs powers of the underlying gauge cou- 
pling g and the Fayet-Iliopoulos coefficient ^. This rescal- 
ing is thus partly responsible for the rise in A^ as a func- 
tion of TV. However, having such large values of A natu- 
rally begs the question as to whether the stability of our 
vacuum towers is in conflict with the assumed pertur- 
bativity of our model. However, as we shall see below, 
there is no conflict between the two, and indeed such 
large values of the rescaled A do not in and of themselves 
undermine the validity of our tree-level calculations. 



D. Perturbativity and Mass Scales 

The results obtained thus far have rested on the as- 
sumption that the physics of our model is perturbative 
at all relevant scales and is therefore accurately approxi- 
mated by tree-level calculations. However, the true vacua 
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of our theory are those field configurations which mini- 
mize the full effective potential V^ff ((/>i), and this includes 
radiative corrections. Thus, the tower of metastable 
vacua which we have presented above is guaranteed to 
be an accurate description of the actual vacuum struc- 
ture of our model if and only if such corrections are small 
and T4fr(0i) ~ Vtrcc, with Vtroc given in Eq. Indeed, 
this must hold within the vicinity of the solution to the 
classical potential for each applicable vacuum index n. 

We will now demonstrate that there is no problem sat- 
isfying all applicable perturbativity constraints, provided 
the gauge coupling g is taken to be sufRciently small. 
Moreover, since small g is also beneficial for stabilizing 
the vacua in the tower, we shall find that there is no 
conflict between the constraints stemming from pertur- 
bativity and those stemming from vacuum stability — 
even as — > oo. 

Note that in this subsection only, we shall revert back 
to our original unrescaled dimensionful parameters in or- 
der to expose the explicit dependence of our physical 
quantities on the gauge coupling g. This will restore fac- 
tors of g and ^ in many of our previous expressions. For 
example, in terms of the original unrescaled dimensionful 
quantities A and ^, our expression for in Eq. p3|) 

takes the form 

^ *2 _ n ^jy) ^ 2 f:2~N /ofi^ 

T{n + y) x[^ + Rn) 

Similar modifications to other expressions occur as well. 

It turns out that a great deal of information about 
the effective potential can be gleaned from non- 
renormalization theorems p^ . For example, in super- 
symmetric field theories, the superpotential is not renor- 
malized by perturbative effects, except via wavefunc- 
tion renormalization. Moreover, this holds true even 
when the superpotential includes non-renormalizable op- 
erators [l^. As a result of such theorems, we expect 
corrections to Vp to arise only at scales near or below 
the supersymmetry-breaking scale g^/£,. However, these 
are typically the energy scales in which we are interested. 
Moreover, corrections to the Kahler potential do not, in 
general, vanish in supersymmetric theories. Thus, it will 
be necessary to discuss both sorts of corrections. 

We begin by addressing radiative corrections to the 
Kahler potential. These arise due to diagrams which 
contribute to wavefunction renormalization of the various 
fields in the theory, and depend on g. For our purposes, 
it will be sufficient to focus on one-loop radiative correc- 
tions to the <f>i propagator; the contributing diagrams are 
then shown in Fig. 1141 Let us begin by considering the 
contribution from diagrams with scalars running in the 
loop. There are N +1 such diagrams, and they result in 
a net contribution 

where rrij is the mass of ipj , where is an arbitrary renor- 




FIG. 14: Diagrams contributing to the wavefunction renor- 
malization of <j)i at one loop. Each of these diagrams is pro- 
portional to . Diagram (a) arises due to the quartic cou- 
plings in the D-term potential, while diagrams (b), (c), and 
(d) arise due to gauge-interaction terms for the scalars and 
the supersymmetrization thereof. 

malization scale, and where 

N 

Tij = QaiQaj ■ (38) 

If Tij were an arbitrary matrix, the expression in Eq. (j37p 
would scale roughly as A^-f 1 for large N, and the theory 
would rapidly become non-perturbative. 

This is not the case, however, due to certain properties 
of Tij which are essentially consequences of the moose 
structure of the model. In particular, all entries along the 
diagonal of this matrix are positive and C(l). Likewise, 
all elements with i ^ j satisfy — 1 < Tij < 0, and the 
sum of elements in any row or column of T vanishes. 
Together, these properties imply that the contribution to 
the renormalization of the Kahler potential from scalar 
loops is essentially independent of N. Moreover, each of 
the additional diagrams in Fig.[T3]yields a contribution to 
the two-point function for (/>i which is roughly of the order 
of Tii (no sum on i). This is also essentially independent 
of N. Consequently, the renormalization of the Kahler 
potential is under control, and radiative corrections of 
this sort can be safely neglected as long as g <C 47r — 
even for very large N. 

The second class of diagrams we must consider are cor- 
rections to the 2iV-field couplings in Vp which come from 
the terms in Vu. In other words, these are corrections to 
the superpotential coupling A which arise from non-zero 
g. The leading such contribution arises at one- loop order 
from diagrams of the sort depicted in Fig. 1151 along with 
additional diagrams in which gauge bosons and gaugi- 
nos run in the loop. Note that similar diagrams were 
examined in Ref. [l7[. In the limit of unbroken super- 
symmetry, of course, these contributions would sum to 
zero. 

Each of the contributing diagrams of the sort pictured 
in Fig. [T51 contains N vertices, and each of these vertices 
contributes a factor of g'^Tij as well as N scalar propaga- 
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FIG. 15: A representative of the class of diagrams which con- 
tribute to the renormalization of the effective 2Af-scalar vertex 
at one loop in the broken phase of the theory, due to quartic 
D-term interactions. The diagram contains A'' vertices (each 
proportional to g^) and A'' scalar propagators (with masses 
proportional to g^^); the overall contribution is thus propor- 
tional to g^. 



tors. In any n-vacuum, the scalar masses are expected to 
be ^ 0{g'^£,R~^), so each diagram is proportional to g"^. 
There are A^! ^ such diagrams, but each is propor- 
tional to TrlTi-^i^Ti^i^ . . .T'i„_iTiv] and hence suppressed 
by a factor of ~ (1/iV)^, where i ^ j. Thus, as was 
the case with the wavefunction-renormahzation calcula- 
tion above, the A^-dependence essentially cancels. Thus, 
as long as g <C 47r, this contribution too can be neglected 
— regardless of the value of N. By the same token, 
contributions to other effective operators which involve 
couplings of various numbers of scalar fields can also be 
safely neglected. 

The final category of radiative corrections we must con- 
sider are corrections to the effective 2Ai"-scalar couplings, 
each of which has the tree-level coefficient A^. Thus, 
these are essentially corrections to the superpotential co- 
efficient A which themselves depend on A. The leading 
contribution arises at one-loop order from diagrams of the 
form shown in Fig. [111 in which 2N — 4 fields are replaced 
by their VEV's, chosen appropriately for a given vac- 
uum. However, along with these contributions we must 
also include the contributions from diagrams of the form 
shown in Fig. [T71 also with VEV's assigned to an appro- 
priate number of external fields. Again, these contribu- 
tions cancel in the limit of unbroken supersymmetry, but 
their contributions can be expected to survive below the 
supersymmetry-breaking scale. 

In the n-vacuum, the one-loop contribution to the 
(Ota^ vertex from a given diagram of the sort 

shown in Fig. [12] is roughly 



167r2 



/(m: 



2 2\ 



n 



1 2 



k=l "k 



N+l^N-n+l"q 



(39) 
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FIG. 16: A representative of the class of diagrams which con- 
tribute to the renormalization of the effective 2Af-scalar vertex 
at one loop in the broken phase of the theory, due to _F-term 
interactions. This diagram represents a 4A'' — 4-scalar ver- 
tex in which 2N — 4 of the scalars on the external lines are 
assigned VEV's, and its contribution is proportional to A^. 
At high energies, when supersymmetry is effectively unbro- 
ken, this diagram is cancelled by the contribution from the 
diagram in Fig. 1171 



'f'l 'f'i+l <^N+1 




^1 *H '^N.l 

FIG. 17: A representative of the class of diagrams which cor- 
rect the effective 2V-scalar vertex at one loop. The net con- 
tribution from this class of diagrams would exactly cancel the 
contribution from the class of diagrams depicted in Fig. I16lif 
supersymmetry were unbroken. 



Here /(m|,m|) is a function of the masses mf and 
whose dependence on these masses is essentially loga- 
rithmic, while q denotes the index of the third scalar 
field (j) whose VEV is missing from the contribution aris- 
ing from the particular diagram in question. Likewise, 
Pjkq is a combinatorial factor representing the number 
of ways of assigning the appropriate VEV's to the exter- 
nal fields. Roughly speaking, we find that Pjkq ~ 2^+^ 
for all {j, k, q). Similarly, there are ~ {N — 2){N — 1) di- 
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agrams which contribute at one-loop order to any given 
2A^-scalar vertex. 

Since the product of VEV's in Eq. ((39|) also appears in 
Eq. p6p . we may rewrite our results directly in terns of 
A^„. Thus, we find that this correction will be small 
compared to the tree-level term in any particular n- 
vacuum as long as 
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■c/ « A^ 



(40) 



Here c/ is an 0(1) coefhcient which embodies the effect 
of including the various f(jnf, m|) functions. 

We recall, however, that we must also satisfy Eq. ^ 
in order to guarantee the stability of our entire vac- 
uum tower. Thus, combining these two results, we find 
that the conditions under which both perturbativity and 
vacuum-stability constraints are simultaneously satisfied 
are given by 



1 < 



A~ 



N 



16^2 



2^+^g'^{N --1){N ~2)cf 



(41) 



Thus, there is no problem satisfying these two inequali- 
ties simultaneously provided 



g' « 



-1){N ~2)cf 



(42) 



We observe that this is also consistent with our previous 
constraint that g <^ At:. 

We conclude, then, that as long as the gauge coupling 
g satisfies Eq. ((42|) and A lies within the range speci- 
fied in Eq. (j4ip . our model will remain perturbative for 
arbitrary N without compromising the stability of the 
vacuum tower. As a result, the tree-level results we have 
presented above are robust against quantum corrections. 
Of course, we observe that the required values of g tend to 
be rather small when N becomes large. However, N need 
not necessarily be taken large for all possible phenomeno- 
logical applications. Moreover, situations in which N is 
taken to be extremely large tend to be higher-dimensional 
deconstruction-type scenarios in which we would natu- 
rally expect our four-dimensional gauge coupling to take 
an extreme value. Indeed, it is not unnatural to expect 
that the fine-tuning inherent in whatever drives N —^ oo 
can also simultaneously drive g — > 0. Unfortunately, the 
details of such a mechanism lie within the full physical 
framework into which such a model is ultimately embed- 
ded, and thus requires a UV completion before they can 
be adequately addressed. 

Our purpose here, however, has been to demonstrate 
that there exists a window in which both perturbativity 
and stability constraints can be simultaneously satisfied 
for any value of A^. As we see from the above discussion, 
this is indeed the case. 

Given these observations, it is interesting to investigate 
the degree to which our model can be considered natural. 
For a given choice of model parameters, and for all A^ > 2, 



this model contains two dimensionful parameters: the 
Wilson-line coefficient A and the Fayet-Iliopoulos term ^. 
Clearly, we can associate a mass scale with each of these 
parameters, defining jix and fi^ such that A = fJ-\~^ and 
^ = /i^. If we assume that ^ and A are generated at the 
same underlying scale /i by the same physics, then our 
model can be considered natural from an effective field 
theory point of view as long as A = c^/i^"^ and ^ = c^/i^, 
where c\ and are both 0(1) coefficients. We shall take 
this to be our definition of naturalness from an effective 
field theory point of view [l^ . 

The question that arises, then, is whether our model 
meets this criterion. Recall that in our model, the partic- 
ular values chosen for A and ^ are constrained by the vac- 
uum stability and perturbativity requirements embod- 
ied in Eq. (|4ip . which in turn depend on the underlying 
model parameters primarily through A^. When written 
in terms of the scales /^a and /ij, this quantity (in rescaled 
variables) is proportional to 
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(43) 



Note, in particular, that this expression contains taken 
to the N/2—1 power. As a result, the extremely large val- 
ues of our rescaled dimensionless A which are required for 
vacuum stability are not in conflict with either perturba- 
tivity constraints or naturalness considerations. Indeed, 
all that is required is that be slightly larger than (but 
still of order) one. 



III. DYNAMICS ON THE VACUUM TOWER: 
CASCADES, COLLAPSES, GREAT WALLS, AND 
FORBIDDEN CITIES 

We now turn to the issue of dynamics within the 
metastable vacuum tower. What will be the pattern of 
tunneling-induced vacuum decays along the entire length 
of this tower? 

Let us begin by recalling the simpler situation that 
arises if we have only two vacua separated by a sin- 
gle saddle-point barrier, with one vacuum state having 
higher vacuum energy than the other. In such a situa- 
tion, the state with higher energy can decay to the state 
with lower energy via instanton transitions, the rate (per 
unit volume) for which may be parametrized as [T^ 
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We will not be particularly concerned with the form of 
the coefficient A. Instead, we will focus our attention 
on the exponent B = Se{4>+t4>-)j usually referred to 
as the bounce action ^ which represents the Euclidean 
action evaluated along the classical path in field space 
which connects , the field-space location of the higher- 
energy vacuum state, to (/<_, the field-space location of 
the lower-energy vacuum state, through the field-space 
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location of the saddle point between them. In general, 
one can evaluate B as in Ref. [20| by approximating the 
potential along the classical path between the two vacua 
as a triangle. In this approximation, the bounce action 
depends on four parameters: A0±, which is the distance 
in field space between the top of the potential barrier and 
the higher-energy (+) or lower-energy (— ) vacuum; and 
AV±, which is the potential difference between the top 
of the barrier and each respective vacuum state. Note 
that in a multi-dimensional field space, the use of these 
results also intrinsically embodies a further approxima- 
tion, namely that the classical path of least action follows 
a trajectory in field space consisting of two straight-line 
segments (one from the higher-energy vacuum to the sad- 
dle point, and the second from the saddle point to the 
lower-energy vacuum). However, this turns out to be a 
reasonably good approximation. 

Calculating B is then relatively straightforward [20| . 
When 



A0_ 



> 



1 



1 



(45) 



with c = (AK-/AV+)(A^+/A^_), the bounce action is 
given by [20j | 
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By contrast, when the inequality in Eq. is not satis- 
fied, the appropriate expression is instead given by (20| 



B = 



9652 \A<j>+^ 
X {-Pi + icPlp. 



+ icPll3+-c^l3l) ,(47) 



where (3± and Rt are given by 
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1/2 



cf3- - P+ 



(48) 



It can be verified that these solutions match smoothly at 
the point where Eq. ([^5]) is saturated. 

Note that the factors of which appear in the denom- 
inators of Eqs. ((47|) and (|48| arise from the fact that we 
are using rescaled energies AV± and field-space distances 
A0-I- in these expressions, in accordance with the discus- 
sion in Sect. II. In the following, we shall take g = 1 
for simplicity in all numerical evaluations of the bounce 
action. 

This is the situation that emerges when there are only 
two vacua to consider. However, in this paper we face a 
situation in which we have a whole tower of metastable 
vacua with many possible pairwise saddle-point solutions. 
The situation we face is therefore significantly more com- 
plicated than that sketched above. 

In order to approach this situation, therefore, we be- 
gin with some preliminary observations. First, we ob- 
serve that if we are interested in transitions between 



an initial vacuum and a final vacuum n/, we need 
only consider the leading quantum-mechanical transition 
amplitude corresponding to the bounce action 

B(ni,nf). Although higher-order quantum-mechanical 
contributions of the form X]n("'/I"') ("l"'^) '^^^ appear 
when there are more than two vacua, such contributions 
are all exponentially suppressed. It is therefore sufficient 
to examine the bounce action B{ni, Uf) itself in order to 
determine the transition rate between two specified vacua 
rii and Uf. 

Second, we observe that in general, a given initial vac- 
uum state rii can decay into all possible final vacuum 
states Uf, where rii < Uf < N — 1. As a quantum- 
mechanical issue, of course, all of these transitions take 
place simultaneously, with rates determined by the corre- 
sponding bounce actions B{ni, Uf). However, once again, 
these transition rates will typically experience huge, ex- 
ponential variations as functions of the possible value n f . 
Indeed, this will be the situation for all < x < 1/2. As 
a result, we shall make a "classical" approximation in 
which each vacuum state rii is assumed to decay to the 
unique final vacuum n/ for which B{ni, rif) is minimized, 
with a rate determined by B(ni,nf). 

Third, in order to evaluate B{ni, n^), we shall need ex- 
plicit expressions for the energies of the rii- and ri/-vacua 
as well as the height of the saddle-point barrier which 
connects them. We shall also require the field-space sep- 
arations between the two vacua and the saddle point. 
While analytical expressions for these vacuum configura- 
tions exist (and were given in Sect. II), we do not have 
analytical expressions for the saddle-point configurations 
except in the "asymptotic" A ^ 00 limit. Therefore, 
although we will not assume that A is actually infinite 
in what follows, we shall assume that A is sufficiently 
large that the asymptotic saddle-point solutions given in 
Sect. II may be utilized without significant error. As we 
have already seen in Sect. IID, this assumption is not nec- 
essarily in conflict with the presumed perturbativity of 
our model; indeed, the approach to the asymptotic limit 
was sketched for the = 3 case in Fig. |4l whereupon 
we observe that the large-A asymptotic behavior emerges 
even for relatively small values of A. The assumption of 
the large-A limit will also have the added advantage of 
removing the free variable A from the subsequent analy- 
sis. 

Finally, we shall deem a metastable vacuum to be "sta- 
ble" if its lifetime exceeds the age of the universe. More 
precisely, we demand a lifetime of such magnitude that 
no decay event would be expected within our Hubble 
volume over a duration equal to the known age of the 
universe « 13.7 Gyr. Using Eq. (|44l) . we can package 
this requirement as a constraint on the bounce action: 



B > 471 + 4 In 
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(49) 



where Afinst = A^^^. When this constraint is satisfied, 
the corresponding vacuum in question is stable on cos- 
mological time scales; when it is not, the vacuum is as- 
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sumed to have decayed sometime in the past. In what 
follows, we will adopt the most conservative assumption 
that Afinst ~ A^pianck, whcrcupon the logarithmic contri- 
bution in Eq. can be ignored. Thus, B « 471 shall 
serve as our critical bounce action for stability. 

Given these assumptions, it is then possible to exam- 
ine the corresponding decay patterns along our entire 
metastable tower. To do this, we need to understand the 
behavior of the bounce actions B{ni,nf) as functions of 
Hi and Uf as we vary our two remaining variables, x ^-^d 
N. As we shall see, there are two principal modes of pos- 
sible behavior patterns ( "collapse" and "cascade" ) which 
will emerge. 



A. Collapse behavior 

To understand how these two different patterns arise, 
let us begin our analysis by revisiting the simple = 4 
case discussed in detail in Sect. IIB. For iV = 4, the vac- 
uum tower contains three minima (the n = 1, n = 2, 
and n = 3 vacua), and hence three different decay tran- 
sitions are possible. The bounce actions B(ni,nf) asso- 
ciated with these three transitions are plotted in Fig. [TH] 
as functions of x- 
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FIG. 18: Bounce actions for the A'' = 4 model, plotted as 
functions of x- In each case we have plotted B{n,n + An) 
where (a) (n. An) = (1,1), (b) (n. An) = (1,2), and (c) 
(n. An) = (2, 1). Also shown (dotted line) is the bounce ac- 
tion corresponding a decay lifetime approximating the age of 
the universe. In general, all bounce actions increase with in- 
creasing Xj ultimately diverging as x ~* 1/2. However, for 
each value of x, we see that the bounce actions corresponding 
to the greatest An are smaller than those corresponding to 
smaller An. 



This figure illustrates several general trends. First, we 
observe the 

• General feature #1: All of our bounce actions van- 
ish as X — > and diverge as x — ^ oo- 

This feature is easy to understand. As x ^ 0, our en- 
tire vacuum tower becomes unstable, whereupon all of 
the possible decays out of any given vacuum state be- 
come essentially instantaneous. Likewise, as x ~* 1/2, 
the vacuum energy differences between any two vacua in 
the tower vanish. There is thus no "driving force" for de- 
cays in the x ^ 1/2 limit, whereupon the lifetime of any 
given metastable state approaches infinity and the states 
become truly stable with respect to instanton-tunnelling 
transitions. 

Second, we observe from Fig. [TH] that B{2, 3) < 
5(1,3) < 5(1,2) for aU x- This implies that the n = 1 
vacuum decays preferentially not to the n — 2 metastable 
state immediately below it, but directly to the n — 3 
ground state. Moreover, if both the n ~ 1 and n = 2 
vacua were somehow initially occupied {e.g., in the dif- 
ferent regions of the universe), we find that the n = 2 
region would decay to the ground state before the n ~ 1 
vacuum region does. 

These observations are examples of two additional gen- 
eral features: 

• General feature #2: A given bounce action 
B{n, n + An) tends to decrease with increasing n if 
An is held fixed. 

• (Nearly) general feature #3: A given bounce ac- 
tion B{n,n + An) tends to decrease with increasing 
An if n is held fixed. (Important exceptions will be 
discussed below.) 

Both of these features are illustrated for the = 20 
model in Fig. [121 where we plot the values of the bounce 
actions B{n, n + An) as functions of n for a variety of 
fixed An. 

These features also have direct physical consequences. 
For example. Feature #3 implies that each metastable 
vacuum in our tower tends to decay directly to the ground 
state of the theory rather than to any other metastable 
state of lower energy. We shall refer to this type of behav- 
ior as a "collapse": each state, one at a time, suddenly 
drops directly to the ground state with n = A^ — 1. As 
we shall discuss below, Feature #3 (and thus the ensu- 
ing collapse behavior) tends to hold for most values of A^ 
and x; indeed, the only exceptions tend to arise in the 
X ^ 0, A^ ^ oo hmit, with both n. An <C A^. Thus, 
collapse behavior tends to dominate along the full length 
of most metastable vacuum towers, and along the lower 
portions of all towers even when x ~* N oo. 

If we imagine situations in which all vacuum states are 
initially populated (e.g., in different regions of the uni- 
verse), the specific collapse pattern along the metastable 
vacuum tower becomes of particular interest. To address 
this issue, we cannot hold n or An fixed; we must vary 
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FIG. 19: Bounce actions for the A*' = 20 theory with x = 1/5. 
For each vacuum 1 < n < 18, we plot B(n,n + An) where 
the different "curves" (from top to bottom) correspond to 
An = 1, 2, 3, 9, 14, 17, 18 respectively. Note that each plot is 
truncated when n + An would exceed rimax = 19. We observe 
that for fixed n, the bounce action B{n,n + An) generally 
decreases with increasing An. As a result, for fixed n, transi- 
tions which maximize An are generally favored. 



400 



2 

m 



200 




100 



FIG. 20: Bounce actions B{n,N ~ 1) with x = 1/5, plotted 
as functions of n for (a) N = 20, (b) iV = 50, and (c) N = 95. 
In each case, the vacua which populate the lower portions 
of each vacuum tower eventually decay by tunneling directly 
to the true ground state. However, the first metastable vac- 
uum to decay in this manner is not the first excited vacuum; 
instead, these decays follow a complicated collapse pattern, 
with different portions of the vacuum tower decaying at dif- 
ferent times. By contrast, for sufficiently large A'^, the vacua 
populating the upper portions of our vacuum towers have life- 
times exceeding the age of the universe, corresponding to the 
critical bounce action B « 471 (dotted line). 



both simultaneously in order to hold n/ = n+An = N—l 
fixed. In other words, we wish to examine the bounce ac- 
tion B{n, N—l) as a function of n. 

This behavior is shown in Fig. [20] for N — 20, 50, and 
95. In each case, we see that the vacua which popu- 
late the lower portions of each vacuum tower tend to de- 
cay first. However, we also see that the first metastable 
vacuum to decay is not the first excited vacuum with 
n = N — 2; instead, these decays follow a complicated 
collapse pattern, with different portions of the vacuum 
tower decaying at different times. For example, in the 
case with A'' = 20 and x = 1/5 shown in Fig. [SOT a). we 
see that the first vacuum to decay into the ground state 
with n = 19 is actually the n = 15 vacuum. The vacua 
then decay sequentially, with decreasing values of n, ex- 
cept that the n = 16 vacuum decays between the n = 15 
and 71 = 14 vacua, the n = 17 vacuum decays between 
the n = 12 and n = 11 vacua, and the n = 18 vacuum 
decays between the n — 6 and n — 5 vacua. 

As TV increases, this behavior persists. Ultimately, 
however, we find that our collapse pattern develops a 
new feature: an upper critical location on the tower 
above which our vacua remain stable on cosmological 
time scales. For example, we see from Fig. [SOf b) that 
for = 50, the collapse pattern begins with the n = 45 
vacuum decaying first. The avalanche of collapses ulti- 
mately spreads up the tower until it reaches the n = 16 
vacuum, at which point it stops. Thus, the vacua with 
n < 16 remain stable. 
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FIG. 21: Bounce actions B{n,N - 1) for TV = 100, plotted 
as functions of n for (a) x ~ 0.1, (b) x = 0.2, (c) x — 0.3, 
and (d) x = 0.33. The critical bounce action (dotted line) 
corresponds to lifetimes exceeding the age of the universe. In 
general, we see that only a narrow band of metastable vacua 
within the vacuum tower will take part in the collapse pattern 
and decay to the ground state; by contrast, vacua above or 
below this band in the vacuum tower will generally remain 
stable. As the kinetic mixing parameter x increases towards 
its maximum value 1/2, this collapse band grows increasingly 
narrow and ultimately disappears. 
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The above results apply for relatively small values of 
X- However, as x increases, our wave of collapses de- 
velops a lower limit as well. This behavior is shown in 
Fig. For example, we see from Fig. l^lTc) that for 
N — 100 and x — 0-3, the avalanche begins with the 
n — 93 vacuum and spreads up the tower to the n = 83 
vacuum before stopping. However, it only spreads down 
the tower to the n = 96 vacuum, where it stops as well. 
Indeed, vacua with n < 82 and n = 97, 98 are stable. In 
general, as the kinetic mixing parameter x increases to- 
wards its maximum value 1/2, this avalanche band grows 
increasingly narrow and ultimately disappears. This is 
consistent with Feature #1 that lifetimes along the vac- 
uum tower diverge as x ^ 1/2. 

We also note from Fig. [20] that 

• General feature #4- For fixed Ui and n/, a given 
bounce action B{ni,nf) tends to increase with N. 

This is a direct consequence of the fact that even when rii 
and n / are held fixed, the corresponding distance in field 
space between the rii- and n /-vacua increases with TV be- 
cause the dimensionality of the field space itself increases 
with N. Since the energies of the n.i- and rt/-vacua are 
iV-independent, the lifetime of the n^-vacuum increases. 



B. Cascade behavior 

As we have seen, it is Feature ^3 which is directly 
responsible for the "collapse" behavior in which each 
metastable vacuum state decays directly into the true 
ground state of the theory rather than into another 
metastable ground state. However, as we shall now dis- 
cuss. Feature #3 is not generally valid, and indeed the 
resulting behavior tends to change rather dramatically in 
the X — !■ 0, — !■ cx) limit. 

To understand the emergence of this qualitatively new 
behavior, let us consider the x — s- 0, iV ^ oo limit analyt- 
ically. We then find that the quantities which determine 
our bounce actions have the leading behavior 
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where we have assumed that rii, rif <C N , and where co, ci 
are C(l) coefficients which generally decrease with rii but 
increase approximately linearly with rif. For example, 
Co = {rif — l)/2 for m = 1. 

The results in Eq. ([50]) are easy to understand. As x ^ 
0, the saddle point between the ni- and n/-vacua shifts to 



join (and thereby destabilize) the n^-vacuum; as a result, 
AV+ and A(/)_|_ both vanish in this limit. However, as 
expected, AV- and A(j)- remain non-zero, even in this 
limit. We also observe from the above results that AV± 
are A^- independent, while all A"-dependence lies within 
A0±. This too is expected. 

Given these results, it is relatively straightforward to 
understand the leading behavior of B{ni,nf) as x ~* 0, 
N ^ oo. We see from Eq. (|50| that there are two limiting 
cases to consider, depending on the value of A^x- For 
Nx/nf <^ Co, the leading term co in Ac/)- dominates. 
In such cases, we find that the inequality in Eq. ([^5]) is 
always satisfied, whereupon the bounce action B{ni, rif) 
takes the approximate analytical form 
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By contrast, for Nx/nf 3> co, it is the Nx/nf term in 
A0_ which dominates. We then find that the inequality 
in Eq. (|45)) is never satisfied, in which case the bounce 
action B{ni,nf) takes the approximate analytical form 
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(52) 

Thus, for any value of A^x? see that the choice of 
whether Eq. ([5T|) or Eq. ([5^ applies depends on how 
A^X compares with con/. However, since co itself grows 
approximately linearly with nf, we see that the first case 
will apply for nf ^ n*r, where 



(53) 



By contrast, the second case will apply for nf n'f- Also 
note that in either case, the actual value of B{ni, n/) de- 
pends on the product A^x- This is not unexpected. Ac- 
cording to Feature #1, the bounce action should vanish 
as X ^ 0, while it should diverge as A^ ^ oo according 
to Feature #4. Thus, when both limits are taken simul- 
taneously, it is not unreasonable that the bounce action 
depend on the product. 

Given these results, we can now consider how B{ni, nf) 
varies with n f for fixed ni in the limit as x ^ Oj A^ ^ c». 
In general, Eq. (j5ip is a rising function of nf, while 
Eq. (j52p is a falling function. As a result, we expect that 
as B{ni,nf) for fixed A^x should generally develop a dip 
(or "trough") centered around n*f as x — ^ 0, A^ — > oo. 
This behavior is shown in Fig. [251 where we plot the 
normalized bounce actions B{\,nf) / B(\,2) as functions 
oi nf in the allowed range 2<n/<Af— 1. In this 
figure, we have taken x = 1/10 and we vary A^ from 
A^ = 10 to A^ = lO'*. As A^ increases, we see from Fig.l^ 
that our bounce-action function indeed develops a deep- 
ening trough whose minimum is approximately located 
ai nf* = \/Nx- This, then, is a counter-example to Fea- 
ture #3. 

This does not, however, eliminate the resulting collapse 
behavior. As long as the minimum bounce action along 
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FIG. 22: Normalized bounce actions _B(l,n/) with x ~ 
1/10, plotted as functions of Uf in the allowed range 2 < 
Hf < N — 1. The curves whose lower portions progress 
from left to right correspond to ordered increasing values of 
N e {1,2,5} X 10*^'^■^^ ending with = 10*. We observe 
that as A'' increases, the bounce actions 5(1, n/) develop a 
deepening "trough" whose minimum is approximately located 
at Uf ~ Uf* = y/Nx- The depth of the trough asymp- 
totes to a finite value as A'^ — » cxa in such a way that the 
minimum bounce action along any such curve corresponds to 
Uf = N -1. 
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FIG. 23: The normalized bounce actions 5(1, n/), plotted 
as functions of n/ in the allowed range 2<n/<7V — 1, with 
iV = 10". The different curves progressing from lowest to 
highest correspond to x = 10"\ 10~^ 10"^, 10~*, and 10~^; 
note that the lowest curve here is the same as the A'' — 10* 
curve shown in Fig. 1221 As x ^ 0, we see that nj 0; thus 
the bounce actions increasingly tend to grow as functions of 
n/, at least for sufficiently small Uf. Also important is the fact 
that for sufficiently small x, the bounce action at n/ = A'^ — 1 
begins to exceed that at n/ =2. For each such curve, the 
minimum bounce action therefore occurs not for Uf — N — 1, 
but for Hf —2. This triggers the onset of cascade (rather 
than collapse) behavior for our vacuum towers. 



any individual curve continues to occur at the maximum 
value Uf — N — 1, the n — I vacuum (and indeed all 
rii-vacua) will continue to preferentially decay directly to 
the Uf = N — 1 ground state. At first glance, one might 
suspect that taking even larger values of N in Fig. [221 
would produce an even deeper trough which would even- 
tually become deeper than the minimum bounce action 
at ny = iV — 1. However, as we see from Fig. [221 the 
depth of the trough actually approaches an asymptote 
as ^ oo, and this asymptotic depth continues to ex- 
ceed the Hf — N — 1 bounce action. Thus our collapse 
behavior remains undisturbed. 

Of course, the above comments pertain to the situation 
in Fig. [221 for which we took x — 1/10. Such a value 
for X — while suitable for illustrating the emergence of 
the trough — is still not yet small enough to alter the 
collapse behavior. The situation changes dramatically, 
however, if we take x even smaller and enter the true 
X ^ regime. For any fixed N, we have seen from the 
above analysis that taking the x ^ limit has the effect 
of sliding 0. We therefore expect that as x ^ Oi 

the initially falling portion of the B{ni,nf) curve should 
disappear, and B{ni, Uf) should actually begin to rise as 
a function of n/ (as long asrif <^ N and N is held fixed). 
This behavior is shown in Fig. 1231 where we have held A^ 
fixed and plotted -6(1, n/) as a function of Uf while we 
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FIG. 24: Same as Fig.[22l except plotted for x = lO"'^ rather 
than X ~ 1/10. We observe that as A'^ increases (approaching 
the true N ^ <x limit), the minimum bounce action begins 
to occur for relatively small values of n/ rather than for rif = 
N — 1. As in Fig. 1231 this triggers the onset of cascade (rather 
than collapse) behavior for our vacuum towers. 
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reduce x from 10~^ to 10~^. Thus, as we see from Fig.[23l 
reducing x has the geometric effect of "uphfting" each of 
the curves in Fig.[22l Indeed, this uphfting is the generic 
behavior that occurs as we enter the true x ~^ hmit. 

Given these results, we see that the net effect of re- 
ducing X is that the behavior shown in Fig. [22] smoothly 
shifts to become the behavior shown in Fig. [M] Indeed, 
these two figures plot exactly the same bounce actions 
5(1, rt/); the only difference between them is that the 
former is calculated with x — 1/10 while the latter is 
calculated with x = 10~^. However, we now see that the 
effect of reducing x has been dramatic: while the smallest 
bounce action i?(l, rif) along any curve in Fig, [^occurs 
for the maximum allowed value nj ~ N — 2, the smallest 
bounce actions B{l,nf) along the curves with large N in 
Fig. l24l now occur for the minimum allowed value Uf — 2. 
Thus, while the top vacuum in our tower prefers to de- 
cay directly to the ground state when x — 1/10, it prefers 
to decay to merely the metastable vacuum immediately 
below it when N is sufficiently large and x = 10^^! 

Since the top vacuum has only decayed to another 
metastable state in the vacuum tower rather than to the 
ground state, the cycle can then repeat. In this case, 
for example, the nature of the subsequent vacuum decay 
now depends on the behavior of B(2,nf) as a function 
of n/ for 3 < ri/ < — 1. However, as long as Ui <C N, 
the same situation persists: the preferred subsequent de- 
cay is into the next-lowest vacuum, and the cycle repeats 
yet again. Thus, what emerges is not a collapse into the 
ground state, but rather a sequential cascade from vac- 
uum to lower vacuum. 

It should be stressed that not all vacuum cascades nec- 
essarily proceed through single- vacuum hops. Indeed, al- 
though this was the result emerging for the curves with 
the largest values of N plotted in Fig. [Ml these curves 
still all correspond to cases with A^x < 1- However, in 
the true x 0, — > oo limit, the size of the cascade 
hops generally depends on the product iVx- This is il- 
lustrated in Fig. [211 where we have plotted the values of 
B{l,nf) for different values of A^X: all while remaining 
completely within the double x ^ 0, A^ ^ oo limit. As 
expected from our above analysis, increasing the value of 
A^X t^as the effect of increasing n*^ and thereby shifting 
the minimum of the curve towards larger values of rij. 
In such cases, the value of n/ which minimizes Bijii, Uf) 
will occur for n/ = rii + An where An > 1. As a re- 
sult, for A^x ^ ^'(1)7 cascades can often proceed through 
larger, multiple- vacuum hops, and indeed non-trivial cas- 
cade trajectories can easily develop. An explicit example 
of such a non-trivial cascade pattern will be presented 
below. 

In general, even in the double x ~* 0; ~^ 00 limit, 
vacuum cascades do not continue all the way down the 
vacuum tower. Instead, beyond a certain point, be- 
comes sufficiently large that the behavior of B{ni,nf) as 
a function of n/ reverts back to the collapse pattern, with 
the minimum bounce action occurring for the ground 
state nf = N — 1. Thus, a system which starts at the 
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FIG. 25: Behavior of 5(1, n/) as a function of of n/ in the al- 
lowed range 2 < nj < plotted for different values of A'^X- 
The curves whose end positions progress from left to right cor- 
respond to ordered increasing values of A^x = IQ''^""'^''^'^'^'^^ , 
where we have held % ~ fixed. As A'^x increases, the 
bounce actions B{l,nf) develop a deepening "trough" whose 
minimum is approximately located at n/ « = V^X- 
a result, the vacuum cascade step size increases as a function 
of iVx- 

top of the tower eventually cascades down the tower to 
a critical level n* at which collapse behavior takes over. 
The subsequent (and final) decay will then be directly to 
the ground state. 

The critical value n* at which cascade behavior be- 
comes collapse behavior generally increases with A^. In 
fact, we have found numerically that 

n* w cA^ (54) 

where c « 2.32 x lO^'^. Remarkably, this relationship 
holds with increasing precision as A^ grows large. More- 
over, this relationship is independent of x- Of course, 
as discussed above, the mere existence of a cascade re- 
gion already presupposes that we are in the x ^ limit. 
However, once x is sufficiently small as to produce cas- 
cade behavior and establish a non-zero value of n*, re- 
ducing X still further will not increase n* beyond this 
value. Indeed, only with a simultaneous increase in A^ 
can an increase in n* be accomplished. 

The result in Eq. implies that for large A^ and 
small X, no more than the uppermost 0.23% of any given 
vacuum tower can experience cascade behavior as opposed 
to collapse behavior. This small value presumably reflects 
the overall numerical coefficients which appear in the ex- 
pressions for our bounce actions in Eq. (|46l) : indeed, we 
observe that c ^ C'(l/327r^) as far as overall scales are 
concerned. Despite its small relative size, however, this 
cascade region of the vacuum tower is of tremendous im- 
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portance for two reasons: 

• For large N, this region gives rise to the vast ma- 
jority of the vacuum transitions that a given sys- 
tem can undergo — indeed this cascade region of 
the tower gives rise to Hterally all possible vacuum 
transitions for our system, except for the final tran- 
sition to the ground state. 

• Likewise, for any N, these vacuum transitions in- 
volve the greatest shifts in vacuum energy that the 
system can experience. Indeed, as we have seen in 
Sect. II, the differences in vacuum energy between 
neighboring vacua grow increasingly small as we 
move down the vacuum tower. Thus, for large N, 
essentially all of the original vacuum energy of our 
system is dissipated through vacuum transitions oc- 
curring in this region of the tower. 

Thus, to summarize: in the double x ^ 0, iV ^ cxi 
limit, we find that a vacuum cascade can emerge near 
the top of our metastable vacuum towers. In general, 
the value of the product Nx governs the size of the hops 
taken in each step of this cascade. Eventually, the cas- 
cade proceeds down the tower until a critical value n* is 
reached, at which point the cascade behavior reverts back 
to collapse behavior. This value of n* increases with N, 
but is essentially independent of x so long as we remain 
in the x ^ limit which was required to produce the 
cascade region in the first place. Moreover, n* is always 
significantly less than N. As a result, there is always 
a substantial collapse region in the lower portion of any 
vacuum tower. 

Finally, of course, we recall Feature #2 which asserts 
that bounce actions which govern the decays of vacua 
near the top of any vacuum tower generally exceed those 
near the bottom. As a result, it is possible that such 
bounce actions near the tops of our metastable vacuum 
towers will exceed the critical value B w 471 required for 
stability on cosmological time scales. As a result, it is 
possible that the uppermost vacua in any vacuum tower 
will be essentially stable. 

C. An explicit example 

To illustrate all of these decay patterns and features 
simultaneously within a single model, let us consider the 
specific example with N = 5000 and x = 2.8 x 10"^. 
Since this choice has x ^ 1 and ^ 1, we see that it 
is precisely such a choice which will yield both a cascade 
region as well as a collapse region. Moreover, since A^x ~ 
0(1), we expect that our vacuum cascade will proceed 
through relatively small steps. 

A plot detailing the vacuum dynamics for this choice 
of parameters is shown in Fig. [211 For each n, we have 
plotted B{n,n') for that value oi n' {n < n' < N — 1) 
which minimizes B{n,n') and which thereby indicates 
the next vacuum along the corresponding cascade trajec- 
tory. For example, we see from this figure that the n = 3 



vacuum decays into the n = 6 vacuum, which in turn 
decays (even more rapidly) into the n = 10 vacuum; this 
in turn decays (even more rapidly) into the n = 15 vac- 
uum, which in turn decays directly into the ground state. 
There are, of course, limits to this cascade region, both 
at the top and at the bottom. For example, the bounce 
actions for the top two vacua exceed our critical bounce 
action B « 471; these vacua, if initially populated, are 
consequently deemed stable on cosmological time scales. 
Likewise, at the n = 11 vacuum and beyond, we enter 
the collapse region in which all subsequent decays auto- 
matically proceed directly to the ground state. 

Nevertheless, within the cascade region between the 
n = 3 and n = 10 vacua, we see that this model contains 
four independent potential cascade trajectories, each of 
which unfolds with increasing speed (i.e., decreasing life- 
times): 
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where 'GS' signifies the ground state. It is therefore only 
an initial condition that determines which trajectory a 
given system ultimately follows. 

Given these results, we can separate the vacua in our 
vacuum tower into three distinct regions on the basis of 
their decay phenomenologies. This is illustrated schemat- 
ically in Fig. [27l At the top of the tower is a stable re- 
gion within which all vacua have lifetimes longer than 
the current age of the universe. This stable region is 
then separated from the remaining unstable regions by 
a "Great Wall" which denotes the border between the 
stable "civilized" world and the remaining "barbarian" 
regions which are affiicted with the omnipresent threat 
of sudden and spontaneous vacuum decay. Moving past 
the Great Wall, the second region is a cascade region 
in which there develops a complex pattern of decays of 
metastable vacua into other metastable vacua. Finally, 
moving further down the tower, the third region is a col- 
lapse region consisting of vacua which decay directly to 
the ground state of the theory. 

Note that this last region may, in turn, be subdivided 
into two distinct (but overlapping) subregions. The first 
consists of vacua (such as the n = {11, 13, 14, 15} vacua in 
the current example) which can be reached at the end of 
one or more decay chains beginning in the cascade region. 
By contrast, the second (which here technically includes 
the n = 3,4,5,7,12 vacua as well as all n > 15 vacua) 
consists of vacua which can never be reached through any 
instanton-tunncUing decay chain. These vacua, which 
may only be populated by some initial condition, collec- 
tively form a "Forbidden City" into which one cannot 
enter from the outside. Indeed, the universe can inhabit 
such a Forbidden City only if it was initially "born" there. 

Thus, we see that constructions of this sort possess 
not only a great number of metastable vacua in addi- 
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FIG. 26: Lifetimes along the difFerent cascade trajectories, calculated for A'^ — 5000 and x = 2.8 x 10~*. For each vacuum n, we 
plot B{n,n') for that value of n' {n < n' < N —1) which minimizes B{n,n') and which thereby indicates the next vacuum along 
the corresponding cascade trajectory. Points are connected according to their sequential trajectories. We also show the Great 
Wall (dotted line) signifying the critical bounce action corresponding to the age of the universe. We see from this plot that 
in this case the top two vacua are stable, while the n — 3 through n = 10 vacua cascade down along four distinct trajectories 
with decreasing lifetimes until they pass the n = 11 vacuum. At this stage, each trajectory decays directly to the ground state 
of our vacuum tower. By contrast, vacua beyond the n = 15 vacuum populate a "Forbidden City": it is impossible to tunnel 
into these vacua from any other locations along the vacuum tower, and the universe can inhabit such a vacuum (for a brief 
time only) if and only if it is initially born there. 
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FIG. 27: A schematic of the vacuum structure of the model shown in Fig. (26] Vacua in the stable region to the left of the 
Great Wall have lifetimes exceeding the age of the universe, while vacua in the cascade region decay to other (lower) metastable 
vacua in the vacuum tower. By contrast, vacua in the collapse region decay directly to the ground state of the vacuum tower. 
Finally, vacua which populate the "Forbidden City" cannot be reached from outside the Forbidden City: such vacua can be 
populated only as an initial condition at the birth of the vacuum configuration. 
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tion to their ground states, but also the potential for a 
highly nontrivial set of vacuum dynamics involving com- 
plicated vacuum cascade/collapse behavior. In this vein, 
it is worth emphasizing that the example illustrated in 
Figs. [26] and [27] represents only one of many possible 
vacuum-cascade scenarios that can be realized in scenar- 
ios of this sort. Furthermore, if we were to relax some of 
the simplifying assumptions inherent in the model out- 
lined in Section [n] — for example, our assumption that 
all nearest-neighbor kinetic-mixing parameters in Eq. ([3]) 
are equal, with all others vanishing — an even wider 
range of possible cascade scenarios would result. A sim- 
ilar possibility exists if A is not taken in the asymptotic 
region: decreasing A tends to decrease the lifetimes of 
our metastable vacua, and thereby enables more rapid 
vacuum transitions to occur. 



IV. SPECTRUM OF THE MODEL 

We now turn to a discussion of the spectra of physical 
particles that arise in the different vacuum states of our 
model. Rather than provide a detailed phenomenological 
analysis of these particles, our main interest is in under- 
standing how their mass spectra evolve as functions of 
the vacuum index n. This will enable us to trace the 
changes in these particle spectra as our system evolves 
down the vacuum tower towards the ground state. 



A. Scalar spectrum 

We begin with the scalar sector, which comprises 
2A^ -|- 2 degrees of freedom: the real and imaginary com- 
ponents of the 4>i. In any given vacuum, the squared 
masses of these states are the eigenvalues of the mass ma- 
trix defined in Eq. (fTT|) . with all fields replaced by their 
expectation values appropriate for that vacuum. How- 
ever, we immediately observe from Eq. (|30p that — 1 
of the (pi receive nonzero VEV's in any vacuum state. As 
a result. A'' — 1 global U{1) symmetries (corresponding to 
global phase rotations of these fields) are spontaneously 
broken, and we expect the spectrum to contain A^ — 1 
Nambu-Goldstone bosons. Thus A^ — 1 of the eigenvalues 
of Eq. ([5D|) will vanish, and the squared masses for the 
remaining A^ -I- 3 scalar degrees of freedom will be pos- 
itive (by definition, since each vacuum is presumed sta- 
ble). Moreover, in any n- vacuum, four of these remaining 
A^ -t- 3 degrees of freedom are the real and imaginary com- 
ponents of the two complex scalars (j>N-n+i and (/)n+i — 
the fields whose VEV's vanish in that vacuum — with 
masses proportional to A^^^~^. The other A^ — 1 degrees 
of freedom are real scalar fields (the real components of 
the (pi for which vf ^ 0), and will have squared masses 
proportional to ^. 



B. Gauge-boson spectrum 

As we have noted above, A^ — 1 of the A^ different 
U{1) gauge-group factors in our model are spontaneously 
broken in each vacuum state by the VEV's of the (pi 
fields. As a result, the A^ — 1 Nambu-Goldstone scalar 
modes discussed above are "eaten" by the gauge bosons 
associated with these broken t/(l)'s, resulting in a gauge- 
boson mass matrix of the form 

N+l 

(A^laugc)afc = .g' ^--Qbr • (56) 

i 

Here we have explicitly written the hats on the charge 
matrices, as in Sect. II, to indicate that they are de- 
fined in the basis in which gauge-kinetic terms take their 
canonical forms. Of course, this matrix (and indeed the 
physical mass spectrum dictated by its eigenvalues) is dif- 
ferent for each vacuum in the tower. However, for all n, 
it has precisely one zero eigenvalue which corresponds to 
the linear combination of gauge bosons associated with 
the single remaining unbroken U{1) gauge group. This 
linear combination, which we will call U{iy with gauge 
field B^, turns out to be a uniform admixture of the 
gauge fields ^7v_„+i, . . . , in the unhatted basis: 

1 ^ 

= -= ^^ (57) 

V a=N-n+l 

Note that this massless gauge boson couples only to the 
complex scalars (pN-n+i and (jjN+i- 

The other A^ — 1 linear combinations of gauge fields 
are nothing but massive Z' vector bosons with masses 
proportional to ^ (and independent of A). Note that the 
mass of each such gauge boson is identical to that of 
one of the A^ — 1 real, massive scalars mentioned above. 
This occurs because all non-vanishing quartic couplings 
among the various (j)i are given (in the unhatted basis) 
by . This results in the massive gauge bosons having 
the same masses as the massive scalars. Indeed, such a 
situation arises as a result of the structure of the D-tevm 
potential in any supersymmetric theory in which there is 
no F-term contribution to the quartic couplings of the 
scalars. 



C. Fermion spectrum 

We now turn to the fermions in our model. These 
consist of the superpartners ipi of the complex scalars (pi 
{i = 1, A^ -f 1), as well as the gauginos Xa associated 
with the U{l)a gauge groups (a = 1, ...,N). As a result, 
the fermion sector of the theory contains a total of 2N+1 
Weyl spinor degrees of freedom. 

Since supersymmetry is broken in each vacuum in the 
tower, one of these spinors must play the role of the Gold- 
stino. Moreover, since supersymmetry is broken solely 
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by _D-terins in each vacuum, we can identify this par- 
ticle with the gaugino superpartner of i?^ in Eq. ([57)1 . 
In supergravity scenarios, this particle is "eaten" by the 
gravitino according to the super-Higgs mechanism. 

All other fermions are physical and acquire Dirac 
masses. In the n vacuum, we find that V'Af-n+i and V'JV+i 
(the superpartners of (pN-n+i and (f>N+i) acquire masses 
proportional to A^^^~^^/^, stemming from superpotential 
contributions. The rest of the fermions acquire masses 
through terms of the form V^ViQaiipiXa -l-h.c. that arise 
from the supersymmctrization of the gauge interactions, 
and the mass eigenstates of this latter group are often 
highly mixed. 



D. Evolution of spectra under vacuum transitions 

Given these results, we now seek to understand how 
these different mass spectra evolve as our system under- 
goes vacuum transitions from one vacuum to the next. 
In order to do this, it will be convenient to separate the 
physical, massive particles in the theory into two distinct 
classes based on whether they do or do not couple to the 
massless gauge field B^^ in Eq. ([57)1 . 

As we have already seen, the particles which couple to 
are the component fields associated with the two chi- 
ral superfields '^N^n+i and \E'Ar+i. These particles have 
masses proportional to and mix only with each 

other. Because these states couple to our only massless 
gauge boson in the theory, we shall refer to such states 
as "couplers". 

By contrast, the remaining states in the theory cou- 
ple only to the extra, massive Z' gauge bosons which are 
orthogonal to B^. As such, they do not couple to the 
massless gauge field itself, and we shall refer to these 
states as "non-couplers". These states are massive and 
highly mixed, with masses and mixings which are rather 
complicated functions of x- Broadly speaking, however, 
these masses fill out a closely-spaced "band" which ap- 
proaches a continuum as ^ oo. Note that the masses 
of the non-couplers are A-independent. 

To help illustrate these features, let us consider the 
spectra of scalars and fermions that arise for the param- 
eter assignments N = 20, x = 1/4, and A = X^q ~ 
1.05 X 10^. These scalar and fermionic spectra are shown 
in Figs. [55] and respectively. In each plot, the masses 
of the "coupler" fields are indicated by the solid curves. 
By contrast, the minimum and maximum masses of the 
remaining fields are indicated by dashed lines, and the 
spectral "band" of masses which they demarcate has been 
shaded. 

These figures highlight several significant features of 
the particle spectra. First, it is evident from these figures 
that the masses of the couplers are highly dependent on 
the choice of vacuum. In particular, it is the lightest 
coupler whose mass vanishes and then becomes tachyonic 
when a given vacuum is destabilized. In this particular 
example, A has been set to the critical value A2o- Since 
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FIG. 28: The spectrum of massive scalars, plotted as a 
function of the vacuum index n for = 20, x = 1/4, and 
A = A20 ~ 1.05 X 10*. (a) The solid lines indicate the masses 
of the two complex scalar "couplers", while (b) the dotted 
lines demarcate the edges of the shaded band within which 
the masses of the remaining scalars lie. Note that the lighter 
coupler actually becomes massless for the n — 6 vacuum; it is 
in this manner that the n — 6 vacuum begins to destabilize 
for A — A20, in accordance with the results shown in Fig. 1121 
for X ~ 1/4. Note that the masses plotted are dimensionless 
rescaled masses, as discussed in Sect. II. 
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FIG. 29: The spectrum of massive fermions, plotted as a 
function of the vacuum index n for A'^ — 20, x = 1/4, and 
A = A20 ~ 1-05 X 10*. (a) The solid line indicates the single 
fermionic "coupler" , while (b) the dotted lines demarcate the 
edges of the the shaded band within which the masses of the 
remaining fermions lie. As in Fig. 1281 the masses plotted are 
dimensionless rescaled masses discussed in Sect. II. 



25 



c 




M 0.5 - 



CD 



FIG. 30: Gauge-boson mass spectra, plotted for the n — 
{1, 5, 10, 15, 19} vacua. For this plot, we have taken A'^ = 20 
and X = 1/4, and the (rescaled) masses are displayed on the 
vertical axis. Note that when n is small, many of the masses 
tend to cluster around a particular value. As n increases, 
however, this value drops; fewer of the masses are clustered 
around this value, and the masses of the remaining gauge 
bosons become more widely spaced. By contrast, the mass 
of the lightest massive gauge boson increases with increasing 
n. Thus, as our system tumbles down the vacuum tower, the 
lightest non-zero gauge boson tends to become increasingly 
heavy. 



X = 1/4, we see from Fig. [T^ that it is the n = 6 vacuum 
which is destabihzed at this value of A, and hence it is 
in the n = 6 vacuum that the hghtest coupler becomes 
massless for this choice of parameters. By contrast, the 
other vacua near the top and bottom of the vacuum tower 
are more comfortably stable for this value of A. As a 
result, the masses of the coupler fields become quite large 
for these vacua. 

The situation is quite different for the non-coupler 
fields. For these fields, we see from Figs. E51 and that 
the boundaries of the non-coupler spectral band remain 
roughly constant as one transitions from vacuum to vac- 
uum. However, it will be noted that both the width of 
the band and the mass of the lightest particle in it in- 
crease slightly when n is near iV — 1 . It is also apparent 
from these diagrams that the lightest massive particles 
in any particular vacuum are generally the lightest non- 
coupler scalar and the lightest massive gauge field, both 
with precisely the same mass. 

Figs. [25] and [25] are calculated for A sitting precisely 
at the critical value A2o- However, it is easy to see what 
happens as we increase A: the solid curves corresponding 
to the coupler masses rescale with A, while the shaded 
bands corresponding to the non-coupler masses remain 
invariant. 



It is important to note that although the maximum and 
minimum of the non-coupler "bands" in Figs. I^andl^ 
remain roughly constant as a function of vacuum index 
n, there nevertheless exists a rich rt-dependcncc for the 
masses of the individual states within the band. This is 
shown in Fig. [30] for the gauge-boson spectrum. When n 
is small, many of the gauge-boson masses tend to clus- 
ter around a particular value. As n increases, however, 
the masses of the remaining gauge bosons becomes more 
widely spaced. By contrast, the mass of the lightest mas- 
sive gauge boson increases with increasing n. Thus, as 
our system tumbles down the vacuum tower, the light- 
est non-zero gauge boson tends to become increasingly 
heavy. 

Given these results, several intriguing phcnomenologi- 
cal possibilities emerge. As we have seen, for any vacuum 
in our metastable vacuum tower, there is only one physi- 
cal, massless field: this is the gauge boson in Eq. (jFT)) , as- 
sociated with the single remaining unbroken U{iy gauge 
group. If this U{iy gauge boson resides in a hidden sec- 
tor — or if it is identified with the photon — there should 
be no difficulties in making this model compatible with 
present experimental observations. In particular, if we 
imagine that all of our broken f7(l)'s are hidden or bro- 
ken at mass scales which exceed accessible energies, then 
only the coupler fields will be readily observable. Since 
the masses of these fields scale linearly with A, we can 
easily adjust them to be unobservable as well. Thus, the 
presence of such large numbers of fields need not neces- 
sarily pose phenomenological difficulties. 

On the other hand, it remains true that the spectrum 
contains a large number of massive gauge bosons, scalars, 
and Dirac fermions whose properties depend on the par- 
ticular vacuum state in question. If such states are ob- 
servable, such a spectrum could lead to a rich and inter- 
esting phenomenology, with a variety of potential impli- 
cations for both collider physics and cosmology. 



V. DEGENERATE VACUA AND BLOCK 
WAVES 

Thus far, we have analyzed our model in the range 
< X < 1/2, and investigated the vacuum towers and 
tunnelling dynamics which result. Indeed, each succes- 
sive vacuum in the tower has a lower energy than the 
previous one, and consequently there exists a net direc- 
tion for dynamical fiow. 

All of this changes, however, if we consider the end- 
point where x — 1/2. Note from Fig. [T]that for any N, 
the point at which x is strictly equal to 1/2 is still within 
our set of allowed values of x- In this case, we find that 
Rn — 2 for all n, and the expressions for the vacuum 
energies Vn and saddle-point energies Vn,n' become inde- 
pendent of their vacuum indices n, n'. As a result, the 
vacuum structure in the asymptotic-A regime discussed 
in Sect. II consists of — 1 degenerate vacua with en- 
ergies (Vn) = 1/2 for all n. These in turn are separated 
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from each other by a set of equivalent, saddle-point po- 
tential barriers of uniform height (F„^„') = 2/3 for all 
n, n' . Furthermore, when x = 1/2, we find that the field- 
space configurations for our vacua and saddle points in 
Eqs. pO]) and ((32)) respectively become essentially identi- 
cal, differing from one another only insofar as a different 
VEV is set to zero in each case. This implies that the 
field-space distance between any pair of vacua is equal to 
that between any other such pair. 

This behavior is shown in Fig. [3lJ where we plot the 
vacua and saddle points of our vacuum towers as a func- 
tion of X- For X < 1/2, we see that our vacuum tower 
"staircase" has a non-zero slope. However, as x ^ 1/2, 
this slope becomes less and less until our vacuum tower 
becomes degenerate, with all vacua and saddle points 
having equal energies and distances in field space. This 
situation is shown in Fig. I3ir c). 
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waves fill out a continuous band. As a result, the vac- 
uum energy of the true ground state of the theory will be 
smaller than the vacuum energy of any individual vac- 
uum. 

The phenomenological implications of such a Bloch- 
wave vacuum structure will be discussed in more detail 
in Ref. [3]. Nevertheless, it is interesting to speculate 
that the vacuum energy of this true ground state might 
actually vanish. If this is possible, we would have a situ- 
ation in which our Bloch-wave vacuum structure actually 
restores supersymmetry. Another possibility is that the 
vacuum energy of this true ground state only approaches 
zero as iV — > oo. In this case, one could potentially obtain 
a ground state with a very small cosmological constant, 
in a manner reminiscent of the proposal in Ref. Q . These 
and other issues will be explored more fully in Ref. Q . 

Needless to say, the existence of such Bloch-wave 
ground states is yet another example of a non-trivial vac- 
uum structure associated with the moose. In some sense, 
the "translational" symmetry that shifts us from one de- 
generate vacuum to the next is nothing but a reflection 
of the underlying translational symmetry of the origi- 
nal moose. We caution, however, that it is not merely 
any moose which gives rise to this sort of band struc- 
ture; it is only one with a precise value for the kinetic 
mixing between nearest-neighbor U{1) gauge factors, as 
well as non-trivial Fayet-Iliopoulos terms. Together, it 
is these ingredients which conspire to produce a vacuum 
structure consisting of an infinite number of degenerate 
vacua, each of which is stable (lacking either tachyonic or 
fiat directions), with a non-zero transition probability be- 
tween them. However, once these features are achieved, 
the ensuing Bloch-wave vacuum structure is inevitable. 



VI. DISCUSSION 



FIG. 31: Vacuum structure in the A'^ = 6 model, plotted for 
(a) X = 0.2, (b) X = 0.4, (c) x = 0.5, and (d) x = 0.54. As 
X increases from zero, we see that the "slope" of our vacuum 
"staircase" decreases, ultimately becoming completely flat at 
X ~ 0.5. For X > 0.5, the vacuum staircase inverts, with the 
lightest states now becoming the heaviest states. 

Taken together, these results imply that that the tran- 
sition rates between all vacua in the "tower" are identical 
for X = 1/2. As a result, all vacuum transitions in this 
setup will occur with identical rates, and it is no longer 
appropriate to employ the classical approximations dis- 
cussed at the beginning of Sect. HI in which we assume 
that only one decay channel dominates the tunneling dy- 
namics. Instead, our system must be treated quantum- 
mechanically. 

The result, however, is clear: the true ground state of 
such a theory is no longer any of the individual n- vacua 
by itself. Instead, what results is an infinite set of Block 
waves across the entire set of degenerate vacua. More- 
over, the vacuum energies associated with such Bloch 



In this paper, we have presented explicit examples of 
two new highly non-trivial vacuum structures that can 
arise in supersymmetric field theories. In particular, we 
have shown that when kinetic mixing among its U{1) 
gauge groups is permitted, an A^-site Abelian moose con- 
struction gives rise to a tower of A^ — 2 metastable vacua, 
each of which involves nonzero VEVs for all but two of the 
"link" fields on the moose, in addition to a stable ground 
state. As A^ is increased, the energies of the existing 
vacua remain unchanged, while new vacua appear with 
smaller and smaller energies, the lowest of which asymp- 
totically approaches zero as A^ — *■ oo. We investigated 
the dynamics associated with transitions between these 
vacua, and found that different regions of the vacuum 
tower could manifest very different instanton-induced 
vacuum-decay patterns. These include "collapse" re- 
gions, in which all vacua decay directly to the ground 
state in a particular order; "cascade" regions, in which 
a potentially complicated decay chain from metastable 
vacuum to metastable vacuum develops; stable regions 
in which the vacuum states have lifetimes exceeding cos- 
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mological time scales; and regions into which instanton- 
induced tunnelling cannot occur. We also explored the 
particle spectra associated with these vacua, and showed 
that in each vacuum state there exists only one physical, 
massless mode — the gauge boson associated with the 
lone unbroken U{1) in that vacuum — along with a tower 
of massive gauge bosons, scalars, and Dirac fermions. 
The massless state is coupled to the rest of the states in 
the tower only by highly-suppressed interactions involv- 
ing extremely heavy particles, and thus can be viewed as 
part of a decoupled hidden sector. 

Needless to say, there are many additional avenues 
through which additional properties of our model might 
be explored. For example, we have already seen in Sect. V 
that the % — + 1/2 limit of our general framework produces 
a Bloch-wave structure for the true ground state of the 
theory. This is another unique vacuum structure which 
has not been explored in the literature thus far, but which 
clearly emerges in models of this type. Models exhibit- 
ing Bloch-wave ground states can be expected to have 
phenomenologies which differ markedly from those based 
on single vacuum states, much as the strong CP problem 
of QCD is a feature uniquely associated with the QCD 
^-vacuum. A detailed examination of the implications of 
such a vacuum structure is forthcoming Q . 

Another salient issue worthy of investigation con- 
cerns the thermal properties of such an infinite tower 
of metastable vacua. In this connection, there are two 
issues which are of paramount importance: that of ini- 
tially populating the landscape, and that governing its 
decay patterns. 

The former issue is critical for determining whether the 
universe tends to start out in the cascade region near the 
top of the tower, in the collapse region near the middle 
of the tower, or at the bottom, at or near the ground 
state. Broadly speaking, vacua which contain more light 
degrees of freedom tend to be preferred in a thermal con- 
text. Occasionally, this means that the universe is far 
more likely to end up in a metastable state than in the 
ground state of a given theory. In the scenarios discussed 
in Ref. 0, for example, it is has been shown [2l[ that 
thermal effects prefer the metastable vacuum. By con- 
trast, in our model, the number of light degrees of free- 
dom in any given vacuum along the tower is essentially 
vacuum- independent. Thus it is not a priori obvious into 
which vacuum state the universe would prefer to settle, 
or whether a (potentially large) number of states would 
emerge as equally likely candidates. This would be an 
interesting area for future research. 

The latter issue is also extremely important, for in 
this paper we have limited our attention to vacuum de- 
cays which proceed through instanton-induced tunnel- 
ing. While this is indeed the whole story at zero tem- 
perature, finite-temperature scenarios contain additional 
sphaleron-\ike processes through which such vacua might 
also decay. This has the potential to significantly modify 
the vacuum dynamics we discussed in Sect. III. For exam- 
ple, while the Forbidden City cannot be entered through 



instanton-tunnelling transitions, it can nevertheless be 
entered through thermal excitations. 

A final question worth exploring concerns the theo- 
retical interpretation of our model. As we have seen in 
Sect. II, our model has as its core an A^-site moose the- 
ory of a sort which is familiar to deconstructionists [23 |. 
This suggests that our model should have a natural five- 
dimensional interpretation in the N —^ oo limit. How- 
ever, our moose theory is complicated by the fact that we 
must introduce non-zero kinetic mixings between nearest- 
neighbor J7(l)'s in order to achieve vacuum stability. One 
wonders, then, about the extent to which this modifies 
our previous five-dimensional interpretation. One possi- 
ble clue comes from the fact that we can "rotate" our 
C/(l) basis in such a way as to eliminate the kinetic mix- 
ing entirely; this occurs at the cost of introducing three 
or more non-zero U{1) charge assignments for each of our 
chiral superfields as well as the introduction of non-zero 
"bulk" Fayet-Iliopoulos terms which are located off the 
moose endpoints. In general, it is known [23] that intro- 
ducing non-uniform gauge charges along the moose is one 
way of realizing warped or other non-trivial geometries. 
It therefore remains to be seen which, if any, warped five- 
dimensional geometry might effectively emerge from our 
model in the N oo limit. This could be important for 
understanding the set of possible UV completions of our 
model. Indeed, some of these completions might also po- 
tentially include gravitational effects, as have been con- 
sidered in other metastable models [H, H^l , 

Clearly, there are also a number of possible applica- 
tions for a vacuum structure of this sort. 

Perhaps the application which most immediately 
springs to mind concerns a potential solution to the 
cosmological-constant problem. Over the past decade, 
several scenarios have been proposed in which a small 
cosmological constant emerges as a consequence of a large 
number of vacua [1, 0, d, @| . Scenarios of this sort tend to 
posit the existence of a "landscape" of vacua with certain 
gross properties, including a vacuum state whose energy 
is nearly vanishing. One then imagines that the universe 
either dynamically tumbles down to this special vacuum 
state, or is somehow born there. 

However, to the best of our knowledge, no explicit 
model with a vacuum structure exhibiting such proper- 
ties has ever been constructed. Moreover, most of the ex- 
isting scenarios have been proposed in the context of the 
string-theory landscape, where the different vacua corre- 
spond not to different minima of the same theory, but 
to separate theories characterized by distinct parameter 
assignments, gauge structures, particle contents, etc. In 
such a context, it is not clear that the instanton methods 
of Ref. [1^ apply. In fact, it is not at all obvious how (or 
even if) transitions between models can occur in such a 
framework. As a result, it is not clear how the landscape 
of vacua can be populated even qualitatively, much less 
quantitatively. 

By contrast, the towers of metastable vacua in our 
model exist at a single point in parameter space — that 
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is, for a single choice of x, and A — and corre- 

spond to different vacuum states within the same theo- 
retical model. The method by which transitions occur 
between one vacuum state of the theory and another is 
therefore well understood, both at zero temperature and 
at finite temperature. As a result, reliable statements can 
be made concerning both the initial, statistical popula- 
tion mechanism for these vacua as well as the dynamics 
associated with transitions between them. 

Regardless of the particular scenario envisioned, it is 
important to note that a solution to the cosmological- 
constant problem requires not only a small energy for 
the vacuum we inhabit, but also that the additional par- 
ticles present in the model be thus far experimentally 
observable. Thus, either the masses of all additional par- 
ticles appearing in that particular vacuum must be heavy 
enough to have thus far avoided detection, or else these 
additional fields must decouple from those of the Stan- 
dard Model and form a hidden sector. Whether the for- 
mer criterion is met in any given vacuum depends sensi- 
tively on the values of our underlying model parameters. 

It is also important to note that we can achieve a true 
ground state with nearly vanishing vacuum energy only 
by fine-tuning the parameter N in our model. Specifi- 
cally, if we wish our ground state to have a very small 
vacuum energy, we will require a large number of sites 
on the moose. This is not surprising, since we cannot ex- 
pect to have a fine-tuned ground state without introduc- 
ing our fine-tuning in some other fashion. Indeed, while 
most previous scenarios merely posit a large number of 
vacua as an initial condition for obtaining a small cos- 
mological constant, we are explicitly realizing this large 
number of vacua as the result of a different, equally large 
number: the number of U{\) gauge factors in an under- 
lying moose. However, we stress that this is a relatively 
small price to pay, since we are obtaining a calculable 
tower of metastable vacua in the process, each of which 
is free of both tachyonic and flat directions, along with a 
true, stable ground state with the small vacuum energy 
we desire. 

There are other potential applications of our model as 
well. For example, one of the major thrusts in recent 
string-theoretic research has been a statistical study [1^ 
of the string landscape [23|. Through such statistical 
studies, one hopes to uncover hidden correlations which 
may be taken as predictions from (or evidence of) an 
underlying string structure at high energy scales. How- 
ever, one important issue that needs to be addressed in 
this context concerns the proper definition of a measure 
across the landscape: in what manner are the different 
string theories to be weighted relative to each other? 
Clearly, the most nai've approach is to count each string 
model equally, interpreting each as contributing a single 
vacuum state to the landscape as a whole. However, it is 
entirely natural that moose theories of the sort we have 
been examining in this paper can appear as the actual 
low-e nerg y (deconstructed) limits of flux compactifica- 
tions [IjI, and as we have seen, such theories give rise 



to infinite towers of metastable vacua. Indeed, for cer- 
tain choices of the underlying parameters in these mod- 
els, literally all of these vacua may be rendered stable 
on cosmological time scales. The question then arises 
as to whether such theories should be weighted accord- 
ing to the infinitely large number of vacua which they 
contribute to the landscape as a whole. In fact, following 
this line further, it becomes crucial to determine whether 
the landscape measure should be defined in terms of dif- 
ferent theories, or in terms of the different vacua they 
contain. Indeed, if the true underlying landscape mea- 
sure is based on vacua, then a theory with infinite towers 
of vacua is likely to dominate any statistical study of the 
string landscape. As such, the phenomenological proper- 
ties of these sorts of models will dominate the properties 
of the landscape as a whole. 

Another potential application of the model described 
here is as a possible hidden sector in a fully developed 
model of supersymmetry breaking. The advantage of this 
would be that a hierarchy between the Planck scale and 
the supersymmetry-breaking scale could arise dynami- 
cally, as a consequence of vacuum tumbling dynamics 
along the tower. In order for this to be viable, how- 
ever, a number of phenomenological issues would need 
to be addressed. For example, while supersymmetry is 
indeed broken in each vacuum state in the tower, it is 
broken only by background values for the I^-terms asso- 
ciated with the various C/(l) gauge groups in the model. 
Consequently, R-symmetry is left unbroken in each vac- 
uum. If the dominant source of supersymmetry-breaking 
is to come from a sector of this sort, then the introduc- 
tion of additional matter will be required in order to ob- 
tain realistic masses for the gaugino superpartners of the 
Standard-Model gauge fields. Likewise, it is possible that 
an i^-term component to supersymmetry breaking could 
be engineered via a a modification of this scenario to in- 
clude additional, vector-like matter in a manner similar 
to that discussed in Ref. |14i |. 

Needless to say, these are only some of the many ap- 
plications such a vacuum structure might have. There 
exist, however, numerous additional possibilities. For 
example, because our construction relies directly on the 
presence of extra C/(l) gauge symmetries, the low-energy 
limit of our setup could have significant implications for 
Z' phenomenology. Indeed, if the coupling between our 
setup and the Standard Model is properly engineered, 
the collider signatures of such a scenario could be rather 
dramatic. Note that preliminary analyses of the mass 
spectra associated with each vacuum in the tower can be 
found in Sect. IV. 

Likewise, this scenario could also have a number of as- 
trophysical and cosmological implications. The vacuum 
cascade we have discussed in Sect. HI involves numer- 
ous first-order phase transitions, and is therefore likely to 
generate topological defects. In particular, given that our 
model contains numerous Abelian gauge factors, there is 
a specific likelihood of generating a network of cosmic 
strings 28]. Moreover, different regions of the universe 
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could potentially exist in different vacuum states along 
the tower, giving rise to domain walls. It therefore be- 
comes an important phenomenological question as to how 
the constraints associated with such domain walls might 
be alleviated. 

In summary, then, it is clear that a number of poten- 
tial extensions and applications exist for the new vacuum 
structures presented in this work. Indeed, as the poet 
William Carlos Williams might well have written, 

So much depends 
upon 

the vacuum 
structure 

glazed with 
excitations 



forming white 
chickens. ll29[] 
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